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Abstract 



This review is aimed to demonstrate the basics and use of for- 
malism of gauge-invariant nonlocal correlators in nonabelian gauge 
theories. Many phenomenologically interesting nonperturbative as- 
pects of gluodynamics and QCD can be described in terms of corre- 
lators of the nonabelian field strength tensors. It is explained how 
the properties of correlator ensemble encode the structure of QCD 
vacuum and determine different nonperturbative observables. It is 
argued that in gluodynamics and QCD the dominant role is played 
by the lowest nontrivial two-point correlator (Gaussian dominance). 
Lattice measurements of field correlators are discussed. Important for 
the formalism theoretical tools, such as nonabelian Stokes theorem, 
background perturbation theory, cluster expansion, as well as phe- 
nomenological applications to the heavy quarkonium dynamics and 
QCD phase transition are reviewed. 
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1 Introduction and motivation 



This review is devoted to the study of gauge theories, specifically gluody- 
namics and QCD 4 . There are now many good books on QCD fl], ||, |], 3, 
which mainly concentrate on perturbation theory. The nonperturbative (NP) 
aspects of QCD (and gauge theories in general) have a somewhat peculiar 
status. On the one hand our very existence is based on one of the most 
important NP phenomenon - confinement. On the other hand perturbative 
QCD disregards confinement and treats quarks and gluons as weakly inter- 
acting. Perturbative QCD is a vast field of investigations which makes most 
of the contents of the QCD textbooks leaving for the NP part of QCD only 
few short qualitative chapters at the end. 

The main unsolved problem of QCD is the quantitative treatment of NP 
phenomena. We still do not have selfconsistent formalism where nonpertur- 
bative quantities, like masses and condensates, are computed systematically 
from first principles, i.e. from the QCD Lagrangian. Nowadays the only reli- 
able way of getting NP results in QCD is numerical simulation of the theory 
on the lattice. This area of research is of great importance, bringing many 
insights into the subject || |J. 

Another source of the NP intuition in QCD is the analysis of quantum field 
theoretical models simplified in one respect or another. These are a huge zoo 
of low-dimensional models (notably QCD 2 ), abelian systems, supersymmetric 
theories and many others. The studies of the corresponding exact solutions 
provide substantial information relevant for physically more interesting cases. 

We have no experimental chance to study the properties of the world 
without light quarks. However, lattice calculations and theoretical analysis 
in the limits where quark degrees of freedom are decoupled show, that the 
property of confinement (i.e. absence of colored particles as asymptotic states 
0) holds in the theory without quarks (pure gluodynamics) as well as in real 
QCD. It suggests to study the physics of confinement starting from quenched 
approximation, as is done, for example, in 1/N C - expansion. The asymptotic 
freedom || |9| and dimensional transmutation phenomena provide the NP 
mass scale, Aqqd and force the theory to leave its perturbative domain at 
large distances. All other dimensionful quantities (such as string tension for 
static external sources, gluon condensate, glueball and hybrid masses and 
widths) in the quenched approximation are proportional to corresponding 
power of Aqcd with some dimensionless coefficients that in principle can be 
computed in NP QCD. 



4 



The chiral effects, which are mostly responsible for physics of lightest 
hadrons, are originated from the NP gluonic sector of the theory as well since 
there is no spontaneous chiral symmetry breaking in perturbative QCD. The 
confinement - chiral symmetry breaking correspondence, however, is very far 
from trivial. We do not touch this set of questions in this review, planning 
to present an extended separate discussion of this subject elsewhere. An 
interested reader can find a lot about the light quark physics in the confining 
vacuum in the gauge-invariant field correlators approach in |L0] . 



An important question is related with possible phases of the theory and, 
in particular, with the NP effects in different phases. Internal parameters of 
the theory (such as number of colours and flavours, light quark masses) and 
external conditions (temperature, chemical potential, number of space-time 
dimensions) both play important role here. They determine the symmetry 
of the ground state at a given point of the phase diagram and correspond- 
ing spectrum of light excitations. Each point of the phase diagram can be 
characterized by the values of condensates, describing the symmetry break- 
ing pattern. The simplest condensates are given by the NP gluon and quark 
condensates: 

Oi — 

(-F^) ; (^> (1.1) 

7T 

The average in ( |1 . 1|) is defined by Q2.5|) below and all perturbative contribu- 
tions are assumed to be subtracted from ([Lip . The NP-generated nonzero 
value of the gluon condensate (T^j corresponds to the breaking of scale in- 
variance at the quantum level (notice, that finite dimensionful parameter 
cannot be generated perturbatively). Nonzero value of the quark condensate 
indicates chiral symmetry breaking M. One can also consider different 
higher order operators as well as nonlocal averages. In the high density - low 
temperature limit the nonzero diquark condensate (ipip) is predicted, leading 
to the phenomenon of colour superconductivity, while chiral symmetry is re- 



stored |T2| , [13], [14]]. The actual values of the condensates also depend on the 
temperature and another parameters of the theory. 

As it has already been mentioned, physics of light quarks will not be 
discussed in the present review. Our main interest is concentrated on a 
systematic method of computing NP effects in the quenched approximation 
from some fundamental NP input, one and the same for all processes. We 
have chosen for this NP input a set of field strength correlators (FC), namely, 

A w.i,.,n.» = Tr {F^ vl {x 1 )^{x 1 x 2 )F^ V2 {x 2 )...F^ nVn {x n )^{x n x 1 )) (1.2) 
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where F^s are the field strength tensors and $(x, y) are the phase factors, 
introduced for the sake of gauge invariance and discussed in details below. 

The systematic use of gauge-invariant variables was adopted for the first 
time by S.Mandelstam in his celebrated paper |Tj|. It was shown, that usual 
Feynman rules of perturbation theory can be recovered from the Dyson- 



Schwinger equations for the averages (|1.2|). The key idea proposed in [16, [L 



18| is to use the gauge-invariant quantities ( |1.2|) as a dynamical input also 
in the nonperturbative domain and to express gauge-invariant observables 
through (|1.2j ) via the cluster expansion. Research activity during the last 10 
years has shown that several QCD processes of phenomenological interest can 
be calculated using the set (|1.2|) as an input. More than that, a systematic 
cluster expansion can be performed, and the first term, corresponding to the 
FC with n = 2 (sometimes called later the Gaussian or bilocal correlator) 
already gives a good qualitative description of most NP phenomena, while 
higher cumulants can be considered as corrections (see Section 2). There 
are evidences that these corrections in the known cases are not large and 
contribute around a few percents of the total effect. Thus one obtains a 
theory with a simple but fundamental input - the Gaussian correlator - and 
the corresponding formalism is called the Gaussian dominance approximation 
or sometimes the Gaussian stochastic model of QCD vacuum. 

The method discussed in the review can be called "fundamental phe- 
nomenology" since it uses correlators (actually the lowest Gaussian one) ( |1.2|) 
as the only dynamical input. The latter is usually given by lattice measure- 
ments. It is worth mentioning that actual profile of the Gaussian correlator 
plays no essential role for the calculation of hadron properties. The necessary 
NP information enters via string tension a, which is an integral characteris- 
tic of Gaussian correlator and thus one has a unique opportunity to define 
the hadron spectrum in terms of one parameter. Note also in this respect, 
that we do not try to reconcile the present approach with any of the partic- 
ular NP QCD vacuum microscopical models, in particular, instanton liquid 
model in this review (see e.g. Jl9[)- As it will be clear from what follows, 



since the quantities ( |TT2] ) are eventually the main input objects, one might 
not in principle be interested what kind of field configurations play the dom- 
inant role for the particular behaviour of correlators. On the other hand, 
most of the dynamical statements of a particularly chosen model about the 
NP content of the QCD vacuum can typically be reformulated as a gauge- 
invariant statements about properties of the correlators ( |1.2| ), especially if 
the model provides some recipee to calculate them nonperturbatively (see 
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also discussion in [ED 



Lattice measurement is a very important source of information on FC. 
First calculation of the Gaussian FC in |21], |22| was followed by more detailed 
investigation in |24], |25|, |26| where also nonzero temperature was explored. 
Both Gaussian and quartic FC have been measured while testing the field 
distribution inside the QCD string [27, 28, 29] and this analysis together with 
independent measurement done in |3(J supports the Gaussian dominance 
picture. 

One should say a few words about the symmetries of QCD vacuum, which 
make it natural to use FC given by (|1.2|). Since the fields as well as Ei and 
Bi are vector-like and Lorentz invariance and color neutrality of the vacuum 
lead to inevitable conclusion that there is no any preferred direction either 
in the coordinate or in the color space in the physical vacuum of QCD, one 
immediately gets 

(F^(x))=0 (1.3) 

This rules out models with constant vector fields. There is an important 
difference with the scalar field condensation in Higgs-like theories. The latter 
is usually characterized by some constant parameter, e.g. value of the scalar 
field in unitary gauge. This is impossible for the vector field condensation 
since the construction of true vacuum state implies an averaging procedure 
over different field configurations in order to get Lorentz and gauge invariant 
ground state. These averaging algorithms, however, may be of different types: 
for example, in models based on ensemble of classical field configurations like 
instanton gas model one performs it by summing over all positions and colour 
orientations of individual instantons, while in the original quantum case the 
summation over properly weighted gauge field configurations is "built in" 
from the beginning in the functional integral approach. 

It is known already for more than 20 years |TTJ that the QCD vacuum is 
filled with strong NP chromoelectric (CE) and chromomagnetic (CM) fields, 
which due to the scale anomaly cause the NP shift of the energy density of 
the vacuum 



16a, 



-(^(0)^(0)) 



Because of the asymptotic freedom (3{a s ) is negative (at least for small a s ) 
thus making the NP shift e energetically favorable. 

Now a crucial question arises, what are the typical scales of the fields, 
contributing to the r.h.s. of (|1.4)? The original proposal made in [jO]] con- 
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sidered these fields as slowly varying in space-time and hence, in the operator 
product expansion (OPE) formalism the condensates enter as constant co- 
efficients in front of the terms with powers of momentum while higher twist 
contributions are suppressed as subleading corrections. 

However, in the course of lattice simulations it was found [[H], [22], that 



typical correlation length T g , characterizing the spatial dependence of gauge- 
invariant bilocal correlator ( |1.2|) is very small, of the order of 0.2-0.3 fm (see 
discussion below). It seems therefore very natural to expect the onset of 
conventional OPE at such large Q 2 that Q 2 » T~ 2 , while at Q 2 ~ T g 2 the 
system crucially feels the nontrivial space-time distribution of the fields and 
usual OPE can break down. 

One can easily see that the method of field correlators (MFC) described 
in the review is a natural development of the QCD sum-rule method, where 
the NP inputs are condensates, i.e. vacuum averages of operator products 
taken at one point. The sum-rule method has natural limits and cannot be 
applied to long distance processes, since information stored in local averages 
is not enough to describe dynamics in the region where stochastic nature 
of the vacuum plays the most important role. Thus higher excited hadron 
states, string-like spectrum and all properties connected to confinement are 
outside of realm of the QCD sum-rule method. 

Another essential ingredient of the method is the idea of Gaussian domi- 
nance. As we shall discuss in Section 3, the Gaussian approximation is a very 
specific model, where all higher correlators ( |1.2j ) are expressible through the 
product of lowest Gaussian ones. Then all correlators with odd n are assumed 
to vanish. A less restrictive approximation is the one where the correlators 
with n = 2, 3 are nonzero and the higher FC factorize into product of lowest 
ones. This approximation is sometimes called the extended Gaussian model. 
It is worth noting that the term "Gaussian" should not be misunderstood. 
The theory under consideration is not a free theory. The n - point Green func- 
tion is assumed to be factorized as a product of bilocal ones, which however 
are not propagators of free fields. The Gaussian dominance has essentially 
dynamical origin, which is not clear yet, but is strongly supported by lat- 
tice calculations and hadron phenomenology. As was mentioned above the 
Gaussian model is quite successful in decribing NP dynamics in all situations 
studied heretofore. Corrections due to higher correlators in the cases studied 
were not large (about a few percent for the static potential). 

Contrary to the QCD sum-rules MFC is applicable both at small and 
large distances. It includes perturbative QCD at small distances with asymp- 
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totic freedom, but also describes the extrapolation to large distances. 

The MFC has been applied to describe the hadronic spectrum for light 
quarks |J3J, heavy quarkonia |35|, [!7|, |39], heavy-light mesons pOfl . 
glueballs and hybrids [O, 03, formfactors and structure functions in the 



proper light-cone dynamics |£| |44j], string formation f45|, diffractive photo- 
production pESf , scattering |J7|, ^8], ^(| |51], Q and also outside of QCD 



to calculate the spectrum of electroweak theory |53fl . 

There exist several reviews which are partly devoted to specific applica- 
tions of MFC, namely, [[54] treats mostly the hadron spectrum, scattering is 
discussed in |55], |56|, |57| , the problem of confinement is reviewed in [jnj, and 
QCD at nonzero temperature and deconfinement in [58], brief overview of 
MFC is presented in ]59" |. 

Another example of applications is the calculation of FC directly in the 
framework of a given field-theory, which is possible in the Abelian Higgs 
model and compact electrodynamics (see detailed review |20| and references 
therein) . 

We do not touch most of these applications in the present review, refer- 
ring the reader to the cited literature. Instead, we concentrate here on the 
general and more fundamental questions, including definition of FC in the 
continuum and on the lattice, the renormalization, connection to OPE, low- 
energy theorems etc. We discuss the basic properties of NP QCD as a testing 
ground for the method, and show how they are described in the method of 
field correlators already in its simplest, Gaussian approximation. 

The plan of the review is the following. Section 2 is devoted to the basic 
formalism of the method: definitions of field correlators, their elementary 
properties, factorization, nonabelian Stokes theorem. The physics of confine- 
ment in the language of field correlators is discussed in details in Section 3: 
cluster expansion, string formation and profile, perturbative properties, low 
energy theorems, confinement of the charges in higher representations of the 
gauge group and some other questions are reviewed. Section 4 treats the 
physics of heavy quarkonia by the method of field correlators. Finally, the 
deconfinement transition and aspects of NP physics at nonzero temperature 
are discussed in Section 5. 
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2 Basic formalism 



In this Section we discuss some basic definitions and general formalism of 
nonlocal gauge-invariant field correlators as applied to SU(N) Yang-Mills 
theory. 

2.1 Definition of field correlators 

For the most part of the review we are working in the Euclidean formalism 
and consider the Euclidean vacuum picture of QCD fields. From the QCD 
sum rules and heavy quarkonia spectrum analysis we know that the gluon 
vacuum is dense, 

G 2 = -^(F^Fp ~ 0.012 GeV 4 . (2.1) 

71 

A dynamical characteristic of such stochastic vacuum is given by a set of 
gauge-invariant v.e.v., which in the nonabelian case have the form 

— — (Tr G^ lUl (xi, x )G fl2U2 (x 2 ,x )...G fl xo)} = A 1)2 ... n (2.2) 

where 

Gfi k v k {xk, xq) = $(xo,Xk)Fn kl , k {x k )<5>{x k ,x ) (2.3) 
and phase factors $ are defined as follows 

$(x,y) = Pexpi / Apdzp (2.4) 
J y 

with the integration going along some curve, connecting the initial and the 
final points. The angular brackets are defined here as vacuum expectation 
value with the usual QCD weight, the exact form of which is given by 

(Q> = J 'dfx(A)e- s ^Q(A) (2.5) 

Here S(A) is the standard Yang-Mills action 

S ( A ) = h\ d ' x Tr F i 
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and the measure dfi(A) includes gauge-fixing and the corresponding Faddeev- 
Popov terms. The fermion determinant is taken to be unity in the quenched 
case considered here. 

The most attractive feature of the nonlocal averages ( |2.2| ) is their gauge- 
invariance, as compared with the case of usual gauge field Green's functions 
(A(x)A(y)..A(z)). There is a serious price to pay, however - - functions 
Ai,2...n depend on the form of the contours, entering Q2.4Q . This dependence 
is to be cancelled in physical quantites, and we will discuss below how such 
cancellation occurs. It will be argued, in particular, that in Gaussian ap- 
proximation contour dependence is a sub leading effect. 

The form of basic quantities (|2.2|) suggests to use the special class of gauge 
conditions, the so called generalized contour gauges , to be discussed 

below. Of particular importance is the case of straight-line contours in ( |2.4|) , 
which corresponds to Fock-Schwinger (F-S) gauge condition [B21 B3]. In this 



case any n-point correlator can be constructed using the invariant tensors 
e^a-p, 5p U , vectors {zi—zq)^ and several scalar functions, depending on relative 
coordinates (zi — z ) 2 , (zi — Zj) 2 , where z is the base point of the F-S gauge, 
which can be placed at the origin. The properties of these functions play a 
crucial role in what follows (see Section 3). 

For the simplest nontrivial 2-point correlator with the points Z\ and Z2 
connected by the straight line phase factors (i.e the point z is chosen on the 
line, connecting Z\ and z 2 ) one parametrizes bilocal correlator conventionally 
as 0,|T7|,|T§: 

tfSw* = y c Ti (( F ^MM^,^)F^ 2 (z 2 Mz 2 , Zl ))) = 

1 ( d d A 

{Z(i2&vii>2 Z U2 S Ul p 2 ) + — \Zt/2^nifi2 Z^^Hi^) I Di(zi Z2)~\- 



2 V<9^i " 2 dz Ul 

This representation does not contain CP-noninvariant part, proportional to 
e Mi^iM2^2 5 which should be included in theories with nonzero #-term. 

Another useful representation comes if one considers correlators of CE 
and CM fields separately : 



^Tr ((E t (x)$(x,y)E J (y)$(y,x))) 
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S ij [D(z) + D 1 (z) + [z i f 



dD x (z) 
dz 2 



+ ZiZ 



dDx(z) 
dz 2 



— Tr ((B i (x)^(x,y)B j (y)^(y,x))) 



S ij [D{z)+D 1 (z) + [zf 



dD x (z) 
dz 2 



dDAz) 



ZiZj dz 2 



where z = x — y. The following notation is useful pl |: 

2 dD 1 (z) 



D\\(z) = D(z) + DAz) + z 



dz 2 



D ± (z) = D(z) + D 1 (z) 



(2.7) 



(2.8) 



The 2-point correlator of the dual field strength F, 



fJ,V 



\^^vafiFafi is con- 
structed from dual functions D, Di in complete analogy with (O), the dual 
functions are related with D and D\ according to 



D(z) = D(z) + 2D 1 (z) + z 



2 dDi(z) 



dz 2 



D 1 (z) = -D 1 (z) 



(2.9) 



It is immediately seen, that on the (anti)-selfdual configurations D\ = 0. 

Let us discuss some basic properties of the functions introduced above. 
First, consider the case of QED p6fl . The gauge field propagator in the 
Feynman gauge takes the form 



(eA^x)eA u (y)) = 5, 



1 e 2 {z) 
An 2 z 2 



(2.10) 



where z = x — y and the square of the invariant charge e 2 (z) (in coordinate 
space) is introduced. Differentiating both sides of (|2.10|) with respect to z 
one finds for {F^ v F pa ) in QED: 



Id (e 2S 



12.11) 



and 



where 



(e 2 F^(x)F^(y)) 



6 1 

^2 ^4 



d(3{e 2 
de 2 



de 2 

z2 % = ^ 2 ) 



(2.12) 
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The formulas ( 2.11 )-( 2.12 ) are obtained in abelian theory. Notice that 
D(z) = to all orders in perturbation theory. The case of Yang-Mills theory 
is different from that mostly in two respects: there are perturbative contri- 
butions to the function D(z) and, which is more important, both functions 
D and D\ acquire NP parts, to be denoted D pert (z) and D np (z) , D 1 ^ (z) , re- 
spectively Pure perturbative contributions to D and D\ are discussed in the 
Section 3, while most of the rest of the review is devoted to the NP parts, 
and we usually omit the subscripts "pert and np. 

The perturbation theory dictates the singular behaviour of D pert } D pert 
at the origin: D pert (z) ~ z~ 4 , if z 2 — > 0. The NP part is normalized to the 
gluon condensate (|2.1|) according to: 

24N C ■ (D np (0) + D7(0)) = (F; v F; v ) = 4vr 2 G 2 (2.13) 

This information is to be used in fitting the lattice data on correlators, which 
we are now in the position to discuss. The first measurements were made 
in quenched SU{2) gauge theory [[H], ^2| . The motivation was to correct 
the predictions of []60], |61j of the shift produced by the gluon condensate on 
the levels of heavy qq bound states: the introduction of a finite correlation 
length strongly reduces the effect on higher levels. The NP part exhibited 
an exponential dependence on the physical distance with a slope T g m 0.16 
Fm. The condensate was also determined with the result G 2 = 0.012GeV . 
Of course the energy scale is conventional, the string tension for SU(2) being 
put equal to the physical string tension cr ps ^GeV 2 . What is actually 
measured are loops made of two plaquettes connected by the phase factors 
along straight lines, with different orientation of the plaquettes. A progress 
was made a few years later when the correlators were measured for 577(3) 
gauge theory in the quenched approximation, i.e. with no dynamical quarks. 
Besides the increase in the computer power available the progress was made 



possible by using cooling technique to study large distance correlations [23 



The configurations produced numerically by use of the correct action are 
filtered by this procedure eliminating short range fluctuations, but leaving at 
the same time intact long range correlations. The principle is the following: 
cooling is done by a local procedure which freezes singular links. It can be 
shown that the maximum distance L which is affected after n steps of cooling 
depends on n in a diffusion-like form : L 2 = D - n. In measuring a correlation 
at distance d a plateau is eventually observed after a number of cooling steps, 
on which the value of the correlator is read, corrected for short distance 
fluctuations. Empirically the procedure requires a minimum distance of 3-4 
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lattice spacings to work. This makes hard to measure correlators at short 
distances. On the one hand the shortest distance observed (say 0.1 Fm) 
should correspond to 3-4 lattice spacings at least. On the other hand to 
avoid dominance of surface effects the lattice must be large enough to cover 



~ 1.5 Fm in distance. In [gg, |26| the correlators was measured for quenched 
SU(3) from 0.1 Fm to 1 Fm in distance. 

At large distances an exponential fit correctly describes the data. At 
short distances a divergent behavior ~ z~ 4 is observed as predicted by per- 
turbation theory. There is no theoretically rigorous way to separate these two 
contributions on the lattice. Different parametrizations were used in [T24[ [26 
to fit the data and they proved to be equally good. One possibility is 



1^) = r ^exp(-|z|/A a ) + A exp(-\z\/T a ) 
a 4 



^D[ at (z 2 ) = ^-exp(-\z\/X a )+A l exp(-\z\/T g ) (2.14) 
The constants A , A\ are related to the gluon condensate G 2 as follows (see 



7l 2 



A + A 1 = — G 2 (2.15) 

lo 

A fit to the lattice data gives for quenched SU(3) (see also Fig.l) 

T s ^0.22Fm ; G 2 = (0.14 ± 0.02)GeV 4 , (2.16) 
i.e. a value large by an order of magnitude than the phenomenological value 



of [11]. Moreover it was found |Ai/A | ~ 0.1. 



In p5[ the measurements were done in full QCD with 4 staggered fermions, 
at different values of quark masses m q = 0.01/a ; m q = 0.02/a. The results 
were in this case 

am q = 0.01 : T g = (0.34 ± 0.02) Fm G 2 = (0.015 ± 0.003) GeV 4 
am q = 0.02 : T g = (0.29 ± 0.02) Fm G 2 = (0.031 ± 0.005) GeV 4 

The value of the condensate can be extrapolated to the physical quark 
masses by the techniques of [|11[] giving G 2 = (0.022 ± 0.006)GeV 4 , an or- 



der of magnitude smaller than the quenched value, and in agreement with 
phenomenological determinations [BO, ED. It is worth mentioned, however, 



14 



that even for full unquenched QCD it is rather small: T g = 0.16 Fm for 
quenched SU(2), 0.22 Fm for quenched 577(3), and 0.34 Fm for full QCD 
with 4 flavours. 

A detailed analysis has also been made of the correlators at finite temper- 
ature |25, 2f|. The 0(4) invariance is lost at nonzero T and five independent 
form-factors exist instead of two at T = 0. 

The main result is that magnetic correlators stay intact in the course of 
the phase transition, while electric correlators have a sharp drop (compare 
Fig. 2 and Fig.3). Above T c the electric longitudinal correlators are very small 
so that the numerical determination is affected by large errors. 

It is of prime importance to stress, that electric condensate, defined by 
the electric parts of fl2.13| ) does not play by itself the role of an order param- 
eter for the confinement-deconfinement transition and, in particular, stays 
nonzero above T c (since NP part of D\ is nonzero). It is only the function D, 
which vanishes at the point of the phase transition and this in its turn indi- 
cates zero string tension (see |3.7|) and hence deconfinement. This behaviour 
of correlators at the deconfinement transition is fully consistent with theo- 
retical expectations. More on the field correlators at nonzero temperature 
and deconfinement transition can be found in section 5. 



2.2 Relations between correlators 



Correlators ( |2.2| ) with different number of F's are not completely indepen- 
dent quantities. The dynamical (Schwinger-Dyson type) equations for them 
are discussed in ||70|| . In this section we are concentrating on the identities 



which must be satisfied due to the Bianchi identities, D^F^ V = 0, which hold 
true for F^ u . The fields are nonabelian and therefore the derivative of the 
correlator with respect to z a contains also the differentiation of the contour 
(see references fFil 



8 ^f ] = i*(z, z')A,(z') + i(z' - z) p /„(z, z'), (2.17) 
where the following notation is used: 

I /ry (z,z')= / da a ^(z^ z + a(z' — z))F fn (z + a(z' — z))- 

J 

■$(z + a{z' - z),z') (2.18) 
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One easily obtains the relation between 2-point and 3-point correlators: 



dD(z) 



"Ml v\a P 



dz 2 



{Tr{F^ Vl {zx)I„(:{zi, z 2 ) (z 2 )^{z 2 ,z 1 )))- 



(Z 2 )l a t(z 2 , Zi)) 



(2.19) 



The formula (|2.19|) is exact and also holds for the perturbative parts of 
the correlators. Let us focus our attention on the NP parts, assuming their 



regular behavior at the origin. One finds [f76 |: 



dD{z) 
dz 2 



z=0 



96N, 



f 



abc I rpa rpb rpc 

\ r H\V\ V\V-2 viVx 



^2.20) 



In a similar way by differentiating any FC one obtains two pieces: "local" 
one, where the field strength operator is differentiated, and "nonlocal" term, 
where the parallel transporters are differentiated. Taking an exterior deriva- 
tive the Bianchi identity eliminates the first term, while the second ones give 
rise to a higher FC. Thus one gets an infinite set of relations between FC. 
For some more complicated relations between FC, e.g. expressing the second 



derivative of D(z) the reader is referred to [81]. 



2.3 Nonabelian Stokes theorem 



The important role in the discussed formalism is played by the nonabelian 
Stokes theorem. There exist many derivations and interpretations of this 
result, see JTi] - |80| . One of the simplest ways |[0| to derive the Stokes 



theorem for the nonabelian fields is to use the so called generalized contour 
gauge 0| (see also [j8~2]), where the vector-potential satisfies the following 
condition t /M (x)A tJj (x) = and is expressed as 



ds gs *p*W,X)) 



(2.21) 



Here the set of contours z„ 



Zfj,(s,x) is such that 



Zfj,(0,x) x „ , z^(l,x) 



Zf,(s,x) = Zfj,(s",z(s',x)) 



(2.22) 



and tn(x) 



ds 
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Consider now the gauge- invariant object — the Wegner- Wilson loop 
34j (which we denote as W-loop throughout this review): 

W(C) = ^Tr Pexp [ij c A^x)dx^j (2.23) 



where the integration goes over the closed contour C. Insertion of (|2.21 ) into 



(|2.23| ) immediately yields the nonabelian Stokes theorem for W(C), 

W(C) = ^Tr Pexp (i J^MG^u, ar )) (2.24) 



where is related to gauge field strength according to (|2.3|) (notice, 
that phase factors along the contours given by z(s,x) equal to unity in the 
gauge (P-21|) ). The surface ordering V in (|2.24|) is induced by linear ordering 
P in (|2.23|) . The reference point Xq is chosen arbitrarily on the surface S 
(for details see |73|] - P0|). Contrary to the abelian case, the surface S in 



(|2.24|) must have the topology of the disk. The inclusion of higher genera 
is nontrivial and requires additional holonomies around fundamental cycles 



The r.h.s. of ( |2.24j ) does not depend neither on the choice of surface S 
nor on the choice of the contours z(s, x) in ( |2.21| ) if the nonabelian Bianchi 
identities D^F^ V = are satisfied. To proceed, one has to consider the 
average of the W-loop {W(C)) over all gauge field configurations. In view 
of fl2.24j ) this average can be expressed through the correlators ( [2.3|) . The 
average (..) and ordering V procedures, which should be applied to W(C), 
are nontrivially related to each other, as we discuss it in the next section. 



2.4 Stochastic vacuum and factorization of field com- 
ponents 

Next we discuss the picture of the Gaussian vacuum. In this model, as it was 
already explained in the introduction, all field correlators are assumed to fac- 
torize into two-point ones. Thus one approximates the complicated ensemble 
of fields of nonperturbative QCD by the one satisfying Gaussian distribution. 
We will discuss in what follows existing lattice evidences for such assump- 
tion. It is very important to mention, that the role of stochastic variables in 
the approach is played by parallel transported field strength variables 
(|2.3j ), not by the gauge field potentials A^. While in perturbation theory one 
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can establish one-to-one correspondence between gauge- covariant and nonco- 
variant descriptions ]15[, Gaussian dominance is gauge-invariant statement, 



and it is not clear, for example, which property of the propagators (A.. A) it 
corresponds to in covariant gauges. 

Whereas in a process with commuting stochastic variables the Gaussian 
factorization is uniquely determined, there are several possibilities for the 
case of non-commuting variables which we discuss now briefly. One way to 
factorize averages (|2.2|) is direct matrix factorization leading to vanishing of 



higher van Kampen cumulants [36 , 37j and direct exponentiation: 



(Tr G(l)...G{2n)) K ~ ((TrG(l)G(2)) (TrG(3)G(4)) . . . (TrG(2n-l)G{2n)) + 

+ all path ordered permutations) (2.25) 

(we remind, that all correlators of the odd power vanish in Gaussian model). 
This scheme is however of not predictive power if we try to evaluate the ex- 
pectation value of two or more W-loops with the help of the basic two-point 
correlator. In order to go further in this direction we have to make asump- 
tions on the factorization of the colour compoments of the fields directly. An 
obvious choice is to assume the following factorization for the components: 

(G ai (l)...G a "{2n)) c sim 

(ja ia2 S a 3 a 4 _ _ _ £an-ia„ ^ G(T)G{2))(Tr G(3)G(4)) . . . (Tr G(2n - l)G(2n)}+ 

+ all path ordered permutations ) (2.26) 

However the difference between the matrix van Kampen factorization 
{.) K and the component factorization (.) c is not fully ordered, and hence it 
would give only a correction to the two-point cumulant, if we had complete 
path ordering. But after applying the nonabelian Stokes theorem the field 
strength matrices are only surface ordered and points in the correlator, which 
are separated by a large number in surface ordering can be physically quite 
close. 

These contributions can destroy the surface law and give rise to a term 
proportional to the 3/2-power of the area. Since surface ordering is much 
less treated in the literature than path ordering we discuss the problem in 
some detail. Let £(£), < t < T be a non-commuting stochastic variable 
describing a centered Gaussian process with correlator: 

(£(*i)£(* 2 )> = (£(* 2 )£(ti)> = c ■ (£ 2 ><K|ii - h\) with 0(0) = 1 (2.27) 
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where c is some numerical constant of the order of 1 and function (j>{x) decays 
with correlation length a. 

In the expansion of the exponential (Pexp[/J"^(t) dt\) we encounter terms 

like 



Pexp[/ T £(t)dt]) = 1+ l T dt x dt 2 (£(tMh))+ 
J o J o Jo 

dh f 1 dt 2 t 2 dt 3 f h dumt,)^))^)^)) + 
Jo Jo Jo 



+ (ati)at3))(at2)au)) + mMuvmMts))) + ■■■ (2.28) 

For the first nontrivial term we obtain: 

f T dh [ h ^ 2 <£(ii)£(* 2 )> = c(e)aT (2.29) 
Jo Jo 

where the value of c depends on the form of the correlation function. The 
next terms crucially depend on the path ordering ( we assume T ^> a): 

f dt x [ h dt 2 f t2 dt 3 f t3 dt 4 (f(*i)f(*2)>(f(t3)f(t4)> = \c 2 (e)\aT) 2 (2.30) 
Jo Jo Jo Jo 2 

f T dh t dt 2 f 2 dh r dU(atiK(t 3 ))(tth)ttU)) = c'(e) 2 a 3 T (2.31) 
Jo Jo Jo Jo 

The term of equation (|2.29|) and the term in equation ( |2.30|) contributes to 

the exponential 

eM\c(0 2 aT] 

whereas the term of equation ( p.30| ) is only of order T and therefore modifies 
the coefficient of T in the exponent by a term of the order c'(£ 2 ) 2 a 3 . The 
details of the factorization can thus always be absorbed in a modification of 
the value of (£ 2 ). 

Let us now consider a stochastic variable on a plane £(u,v) with simi- 
lar properties as the one discussed above (|2.27| ) and investigate the surface 
ordered expectation value: 



(V exp / du dv £(w, v)), 
Jo Jo 



where V stands for surface ordering analogously to (|2.24|) . The first non- 
trivial term in the expansion of the exponential, analogous to fl2.29|) yields a 
contribution proportional to 

{i 2 )a 2 RT 
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independent of surface ordering. But the terms containig the expectation 
value of four values of the stochastic variable can behave quite differently: 
Let again fx, / 2 , /j, fi denote four ordered points on the surface. Then the 
integral over the ordered expression 

yields (£ 2 ) 2 a 4 (.RT) 2 which together with the former term exponentiates and 
contributes to the area law 

ex P [i C (£ 2 )a 2 ^n 

but the ordering 

is not simply a correction to the first term: if fx and fa are close but in 
different rows the points f 2 , f± can vary over the full row, thus yielding a 
factor of the order T , and the final contribution will be of the order 

(£ 2 } 2 a 3 T(RT) 

and thus not fit into an exponentiation. So in contrast to the path ordering 
where different clusterings, as long as they agreed in the leading (ordered) 
term exponentiated to an perimeter (area) law, in surface ordering we need 
a fine tuning in order to get the latter. 

This is a reflection of the rather artificial feature of the Gaussian model 
that we have always to treat correlators where the reference point is fixed 
from the beginning and hence might be quite distant from the connection 
line of the two external points. 

We can exclude this unphysical contributions if we specify the factoriza- 
tion prescription ( |2.26|) by the following ordering rule f45fl : if the expectation 
value is expanded as a power series in the two-point correlators only the lead- 
ing powers of the ratio of the correlation length and the linear dimensions of 
the loops have to be taken into account in each order. 

It is easy to see that the factorization ( 2.26 ) together with this rule leads 
in the case of one loop to the same area law as the matrix factorization (|2.25 ) . 
It will be referred to in the following as the modified component factorization 
in the model of the stochastic vacuum. 



It is instructive to mention here a proposal of [58] based on a domain 
picture of the QCD vacuum in which these problems do not occur and the 
formal modification of the component factorization is justified from a physical 
point of view. 
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3 Confinement and field correlators 



In this chapter we discuss the following question: how confinement is related 
with the properties of FC, what is the physical mechanism behind it and how 
the cluster expansion operates in the quantitative description of confinement. 
In what follows the role of the order parameter of confinement is played by 



the averaged W-loop [83, EM (we assume quenched approximation) 



(W(C)) = ^(Tr Pexp(z f^dxj) = exp[-V(R)T] (3.1) 

This is the true order parameter only in case when light quarks are absent 
and V(R) is the interquark potential for (infinitely) heavy static quarks. 
Confinement for light quarks and gluons is much more intricate phenomenon, 
which is, however, ensured by the same FC as for heavy quarks (see ||34|| ). 

This Section has the following structure: confinement in terms of the clus- 
ter expansion is explained in subsection 3.1, corresponding physical aspects 
and confinement of the charges in higher representations are reviewed in 3.2, 
abelian projection method in terms of field correlators briefly commented in 
subsection 3.3, perturbative structure of the cluster expansion is explained 
in subsection 3.4, the string profile is discussed in subsection 3.5, and, finally, 
subsection 3.6 is devoted to the explanation of the model consistency with 
the low-energy theorems and related questions. 

3.1 Confinement and cluster expansion 



The nonabelian Stokes theorem and cluster expansion theorem applied 
to the W-loop yields the following result: 

(W(C)) = -^(tt Pexp(z£^(x)^) 
= jjr (ti V exp(i J^G^u, x )da tll/ (u)) 

= — Tr exp - J da{l)..da{n)((G{l)..G{n)))\ (3.2) 

where we have suppressed the indices, da(k)G(k) = G^iuk, xo)da^ u (uk), and 
we have used irreducible cumulants. Due to colour neutrality of the vacuum 
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all cumulants are proportional to the unit matrix in the colour space which 
makes unnecessary colour ordering in the r.h.s. of ( |3.2| ). 

The irreducible correlators ((••)) are straightforwardly extracted from 
( |3.2j ) comparing terms of the same order in G. Notice, that there is only 
one ordering in the expression (|3.2j ) and this ordering is induced by the P- 
ordering in the original loop (|3.1| ). 

The most important for us is the bilocal correlator, which is defined as 

\G f j, ll > 1 (xi)G l j i2V2 (x2)) = \(G nwixijG \i2v 2 { x 2))) + (G ^^(xi)) (G ^ 2U2 (x 2 )) 
Since (G(k)) = 0, one can rewrite ( |3.2| ) identically as 

(W(C)) = — Tr exp ( ~J J da^v^da^^K^p^u.v, S)\ (3.3) 
c V s s ) 

where we have defined the global correlator A(u,v, S) [|39| 

-^A FrfW (u,v, S) = -~ ((G flv (u,x )G pa (v,x )))+ 

+ da{3)..da{n) ({G^{u,x )G pa {v,x )G(3)..G{n)))+ 

n=3 n - J 



+ perm(l, 2, ..n) 



(3.4) 



and perm(l, 2, ..n) stands for the sum of terms with different ordering of 
G(u,x ) = G(l) and G(v,x Q ) = G{2) with respect to all other factors G{k). 
Since W{C) does not depend on the shape of the surface S, the dependence 
of the global correlator on its arguments is such, that r.h.s. of ( |3.3|) is inde- 
pendent of the choice of 5* but depends on the contour C. This circumstance 
explains the termin "global correlator" used for A(u,v, S) in contrast to lo- 
cal correlators which enter in the r.h.s. of (|3.4|) (note however, the specific 
meaning of the definition "local" - correlators ((G(l)..G(k))) depend on the 
points zi, .., z k as well as on the paths, entering via transporters $(z, x )). 

If one integrates over all but one variables in (^.3[), it leads to the effective 
integrand 

Q^ v (u,S) = - J da p(7 (v)A^ pa (u,v,S) (3.5) 
s 

In the confining phase one expects for large contours C the minimal area law 
of W-loop, which implies that Q pu in this limit does not depend on the point 
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u when S is the minimal area surface and simply coincides with the string 
tension cr, while for the arbitrary surface one can identify Q^ u as 

Qtiv{ u i S) = P^v • o (3.6) 

where P^ u projects onto the minimal surface. Conversely if Q^ u does not 
have constant limit for large S then the area law of W-loop does not hold. 

To calculate Q^iu, S) one can for simplicity take xo in ( |3.4| ) to coincide 
with u. Then the lowest order FC in ( |3.4|) depends only on two points (as in 
( |2.6j )). In what follows we concentrate on contributions to cr and therefore 
take for simplicity a planar contour C with the minimal surface S lying in 
the plane. 

Insertion of (|2.6|) into (|3.4j) yields the area law of W-loop with the string 



tension cr 

(W(C)) = exp(-aS miri >. 



a W = ly d 2 xD{x) (3.7) 

Note that D\ does not enter cr, but gives rise to the perimeter term. It 
can be clearly seen from the definition ( |2.6j ) since D\ by definition contains 
derivatives. One concludes therefore, that the presence of nontrivial non- 
perturbative contribution to the function D(x) such, that Q3.7p is nonzero 
but finite implies confinement in Gaussian approximation. The lattice data 
described above support this since D(x) demonstrates exponential falloff at 
large distances and hence 

a ~ G 2 • T g 

Till now we have discussed confinement in terms of the lowest cumulant 
- D(x), which is justified when stochasticity condition ( |3.9| ) is fulfilled and 
D(x) gives a dominant contribution. Let us now turn to other terms in 
the cluster expansion ( |3.2|) . It is clear that the general structure of higher 
cumulants is much more complicated than ( |2.6| ) , but there always Kronecker- 
type terms D(xi, x 2 , x n ) Jl <5 WMfe present similar to D(x — y) in ( |2.6D and 
other terms containing derivatives like Di(x — y). The term D(x\, x 2 , x n ) 
contributes to string tension. This means that the string tension in the 
general case is a sum, 

oo „ 

a = ^ (n) , ^ (n) ~ / ((G(l)G(2)...G(n)))da(2)...da(n) (3.8) 
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If one assumes, that all gluon (irreducible) condensates and the corre- 
sponding correlation lengths have the same order of magnitude, then the 
expansion in ( |3.8| ) is in powers of (Xy^G^T 2 ) , and if this parameter is small, 



\^G 2 T 2 g < 1, (3.9) 
one gets the limit of Gaussian stochastic ensemble where the lowest (quadratic 



in G^y) cumulant is dominant. The value of VG2 in ( |3.9| ) is a typical 
scale of vacuum fields, which can be estimated from the gluon condensate, 
\[G~2 ~ (0.6 GeV) 2 , while A ~ 1 takes into account possible cancellations in 
the magnitude of the cumulant, and in this way characterizes the deviation 
of the vacuum ensemble from Gaussian one. Taking into account, that \[G~2 
in W(C) enters on a plane with fixed (ij), i.e. \fG~2 ~ J (G 2 j), while the 
gluonic condensate is a sum J2ij(G 2 j), an d the numerical value of T g from 
lattice calculations T g ~ 0.2 — 0.3 fm, one obtains as the parameter of cluster 
expansion \G2Tg ~ A • 0.1. 

It is of prime importance to stress, that the relative contribution of dif- 
ferent terms in the r.h.s. of the expression Q3.8|) depends on the chosen 
integration surface. Since the l.h.s. of ( |3.8| ) does not depend on the surface, 
the lowest cumulant dominance is expected to hold for the minimal surface 
in Gaussian dominance scenario, while contributions from higher terms are 
to become important for arbitary nonminimal surface in order to cancel ar- 
tificial surface dependence. In can be said therefore, that the actual value 
of the string tension is determined by the total sum ( |3.8| ) in any case: the 
stochastic one (actually governed by Gaussian average for the minimal sur- 
face) and coherent one (when all terms are of the same order even for the 
minimal surface). 

Summarizing, when the stochasticity condition ( |3.9j ) holds, the lowest 
term, a^ 2 \ dominates in the sum; in the general case all terms in the sum 
( |3.8j ) might be important. Therefore the smaller the correlation length T g 
is, the better Gaussian model of confinement works. It is worth mentioning, 
that this limit is opposite to the sum rule limit (sometimes called Leutwyler- 
Voloshin limit), where the condensate is assumed to be constant in the po- 
sition space, i.e. using the language of MFC T g — > 00. From this point of 
view the analogy between Gaussian dominance in the stochastic picture and 
"vacuum insertions" procedure in OPE is not straightforward. 
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3.2 Physical mechanism of confinement 

The described formalism is valid both for abelian and nonabelian theories, 
and it is interesting whether the area law and nonzero string tension could 
be valid also for QED (or U(l) in the lattice version [|90|). To check it let us 
apply the operator |e M i, aj 3gf- to both sides of (|2.6|) f76fl . 



In the abelian case, when $ cancel in ( |2.2| ) and F = G, one obtains 

d a {F a(3 (x)F X a(y)) = e X a 1 pd 1 D{x - y), (3.10) 

where F a/3 = ^e a/ 3^ u F^ u . 

If magnetic monopoles are present in abelian theory (e.g. external Dirac 
monopoles) with the current j^, one has 

d a F a/3 (x) = j p {x) (3.11) 

In absence of magnetic monopoles (for pure QED) the abelian Bianchi iden- 
tity requires that 

d a F af3 = (3.12) 

Thus for QED (without magnetic monopoles) the function D(x) vanishes due 
to ( |3.11| ) and hence confinement is absent (compare with ( |2 . 1 1 ) . 



In the compact version of U(l) magnetic monopoles are present and the 
lattice formulation of the method would predict the confinement regime with 
nonzero string tension ||91|| . The latter can be connected through D(x) to the 
correlator of magnetic monopole currents. Indeed, multiplying both sides of 
( |3.10| ) with ^Xa-yS-^- one obtains 



(Mx)j s (y)) =(^~- ^-6 ps - D(x - y) (3.13) 



dx a dy a dxp dy 5 



This identically satisfies monopole current conservation: applying or 



to both sides of ( 3.13 ) gives zero. More on MFC in abelian theories can be 
found in [M. 

Let us turn now to the nonabelian case, again assuming stochasticity 
condition (|3.9| ). Applying as in the Abelian case the operator \e^ ua f3-£- to 



the r.h.s. of (|2.6j) , one obtains the relation given above in (|2.19|) Especially 



simple form of ( p. 19 ) occurs when tending x to y; one obtains equation ( p. 20 ) 



which can be rewritten as 



dD(z) 



dz 2 



{EtE b Ml)e i3k f abc (3.14) 
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Thus confinement (nonzero a due to nonzero D) is related in nonabelian case 
with the presence of purely nonabelian correlator. 

To see the physical meaning of this correlator one can visualize magnetic 
and electric field strength lines (FSL) in the space. Each magnetic monopole 
is a source of FSL, whether it is a real object (classical solution or external 
object like Dirac monopole) or lattice artefact. In nonabelian theory these 
lines may form branches and e.g. electric FSL may emit a magnetic FSL at 
some point, playing the role of magnetic monopole at this point. This is what 
exactly nonzero triple correlator (Tr FFF) implies. Note that this could be a 
purely quantum effect and no real magnetic monopoles (as classical solutions 
or lattice artefacts) are necessary for this mechanism of confinement. 

The formula ( |3.14| ) demonstrates, that in strict Gaussian approximation 
D(z 2 )' = 0, since the r.h.s. is proportional to the condensate of the third 
power. The extended Gaussian approximation, which keeps only 2- and 
3-point correlators allows to have nonzero D(z 2 )' = 0. It does not auto- 
matically mean, however, that the contribution of the triple (and higher odd 
order) condensates themselves to the physical quantities (like string tension) 
is comparable with that of 2-point one (see (|3.8| ) and discussion below (|3~9|) ). 
It should also be noticed, that despite the first derivative of the function D(z) 
vanishes at the origin in Gaussian vacuum, the second and higher ones of even 
order do not (see [jSTJ]) . 

We conclude this section with discussion of confinement for charges in 
higher representations. Most of the models of NP QCD are designed to 
describe confinement of colour charge and anticharge in the fundamental 
representation of the gauge group SU(3), i.e. the area law for the simplest 
W-loop and hence linear potential between static quark and antiquark. At 
the same time Gaussian dominance scenario naturally describes lattice data 
on the interaction between static charges in higher SU(3) representations. 

We define static potential between sources at the distance R in the given 
representation D as: 

V D {R) = - lim i log(VMC)), (3.15) 

1 — >oo ]_ 

where the W-loop Wd(C) for the rectangular contour C = R x T in the "34" 
plane is given by 

(W D (C)) = (ti d Pexp ^ J A a ^T a dz)j ^ (3.16) 
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The SU(3) representations D = 3, 8, 6, 15a, 10, 27, 24, 15s are character- 
ized by 3 2 — 1 = 8 hermitian generators T a which satisfy the commutation 
relations [T a ,T b ] = if abc T c . One of the main characteristics of the represen- 
tation is an eigenvalue of quadratic Casimir operator C D , which is defined 
according to C$ = T a T a = Cp ■ 1. Following the notations from W2 



we 



introduce the Casimir ratio dp = Cp/Cp, where the fundamental Casimir 
Cf = {N 2 — 1)/2N C equals to 4/3 for SU(3). The invariant trace is given by 
Tr D l = 1. 

The potential (^.15|) with the definition (|3.1(j|) admits the following de- 



composition 

V D (R) =d D V {2 \R)+d 2 D V i ' 4 \R) + ..., (3.17) 

where the part denoted by dots contains terms, proportional to the higher 
powers of the quadratic Casimir as well as to higher Casimirs. It is worth 
mentioned, that the term (R) is a sum of the bilocal correlator contribu- 
tion and linear Casimir contributions coming from higher correlators. 

The fundamental static potential contains perturbative Coulomb part, 
confining linear and constant terms 

V D (R)=a D R-v D -^§ (3.18) 

K 



The Coulomb part is now known up to two loops |)5|, |96| and is proportional 



to C D . The "Casimir scaling hypothesis" @, see also H ||, [10(| |10T 
declares, that the confinement potential is also proportional to the first power 
of the quadratic Casimir Cp, i.e. all terms in the r.h.s. of (|3.17| ) are much 
smaller than the first one. In particular, for the string tensions one should 
get a D /a F = d D . 



This scaling law is in perfect agreement with the results found in 



sec 



also ||102|| ). It is easy to see that Gaussian cumulant in ( |3.16| ) is expressed 



through Cd and representation-independent averages as 

Tr D (F(1)F(2)} = (F a (l)F a (2)} = ^ <F*(1)F«(2)>, (3.19) 

so Gaussian approximation satisfies " Casimir scaling law" exactly. This fact 
does not depend on the actual profile of the potential. It could happen, 
that the linear potential observed in |92| is just some kind of intermediate 



distance characteristics and changes the profile at larger R (as it actually 
should happen due to the screening of the static sources by dynamical gluons 
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from the vacuum) The coordinate dependence of the potential is not directly 
related to the Casimir scaling, and can be analized at the distances which 
are small enough not to be affected by the screening effects (see discussion 
1041 , MPST 



m 



The quantitative analysis ||f05| of the data from |92| shows, that scal- 
ing holds with an accuracy of a few percents. The Fig. 4, taken from the 
paper |92| shows the potential between static sources in the different rep- 
resentations of the gauge group SU(3). Lines on this figure correspond to 
the fit of the fundamental static potential multiplied by the corresponding 
ratio of the Casimir factors. The picture demonstrates nice agreement with 
the Casimir scaling hypothesis. It is worth mentioning that most of the 
microscopic models used in NP QCD (e.g. instanton model, center vortex 
model, abelian dominance hypothesis) encounter difficulties in explanation 
of Casimir scaling [ |104| , |106| . 



In the MVC the Casimir scaling has two important features. First, it is 
the direct consequence of the Gaussian dominance hypothesis since bilocal 
correlator provides the exact Casimir scaling. On the other hand, it implies 
the cancellations of C 2 D - proportional terms and higher ones in the cluster 
expansion of ( |3.16| ). Physically, it means the picture of the vacuum, made of 
relatively small colour dipoles with weak interactions between them. One can 
imagine two possible scenario. According to the first one, Casimir scaling is 
just a consequence of Gaussian dominance. It this case any higher cumulant 
contributes to physical quantities much less than the Gaussian one due to 
dynamical reasons. There is also the second possibility, when each higher 
term in the expansion fl3.16|) is not small, but their sum demonstrates strong 
cancellations of Casimir scaling violating terms. These pictures are in close 
correspondence to the stochastic versus coherent vacuum scenario, discussed 
above. The clarification of the situation in future lattice measurements is 
important to understand the NP structure of QCD vacuum in gauge-invariant 
terms. 



3.3 Perturbative analysis of correlators in gluodynam- 
ics 

Perturbative expansion of FC has specific features due to the presence of 
parallel transporters in the definition of FC, which should also be expanded 
and generate new class of diagrams. In this respect FC are similar to W- 
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loops. At the same time since FC contains a product of operators one expects 
at small distances singular contributions as in the OPE coefficient functions. 
Let us briefly describe the renormalization procedure for W-loops. It was 
studied both for smooth and cusped contours ||112j , |113|| . In the first case the 
result is: 

W{C) = Z ■ W ren (C) (3.20) 

where the (divergent) Z-factor contains linear divergencies arising from the 
integrations over the contour 

Z ~ exp I c J dx^ J dy u D^(x - y) J ~ exp (c — \ (3.21) 
V c c I ' ' 

where L is the length of the contour, a - ultraviolet cutoff and c - numer- 
ical constant not important for us here. All logarithmic divergencies are 
absorbed into the renormalized charge g ren {n) defined at the corresponding 
dynamical scale \i. It is worth noting that a single diagram in the pertur- 
bative expansion of the r.h.s. of Q3.20Q typically gives contributions to both 
Z-factor and W ren (C), i.e. contain linearly and logarithmically divergent 
parts. Hence the factorization in the r.h.s. of (|3.20|) has double-faced nature 



- it is factorization of the diagrams together with the separation of contri- 
butions, coming from dangerous integration regions (i.e. linearly divergent 
terms) for each diagram itself. As a result contour divergencies can be sep- 
arated as a common Z-factor, which physically renormalizes the bare mass 
of the test particle, moving along the contour, which forms W-loop. There- 
fore the W-loop can be calculated or measured on the lattice only modulo 
divergent Z-factor. However if one extracts the physical quantities from the 
W-loop average (e.g. by forming the Creutz ratio) then only logarithmic 
charge renormalization causes an observable effect. 

If the contour possesses cusps (and selfintersections) situation becomes 
more complicated. Namely, each cusp leads to its own local Z - factor, 
depending on the cusp angle (but not on other characteristics of the contour) . 
The expression ( |3.20| ) for the loop with k cusps having angles 73 is to be 



modified in the following way ||1 12[ |113 



W(C) = l[Z l ( lt )-W ren (C) (3.22) 
i=i 

The mixing between linear and logarithmic divergencies, mentioned above 
has some specific features for the loops with cusps. Namely, each cusp in 
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addition to the linearly divergent factors, as in the smooth contour case 
introduces logarithmically divergent terms which have nothing to do with 
the coupling constant renormalization and can be subtracted by their own 
Z(ji) -factors. There are two types of terms of that kind, those, depending 
on the cusp angle and others which do not depend on it. The latter play 
important role in the perturbative expansion of FC. To see this, note, that FC 
can be easily obtained from the W-loop functional. Indeed, differentiating 
( |3.2j ) with respect to 5a flll (xi) and assuming that the contour C connects the 
points u, v with the reference point xq along the straight lines, one finds for 
the bilocal correlator: 

1 5 2 (W(C)) 
-Tr ({G, v (u,x )G pa (v,x ))) = -=- ) \ }) - (3.23) 

and analogous expressions for higher cumulants. Note, that this contour 
generally has at least two cusps (typically even four, if the reference point x 
does not coincide with one of the correlator arguments). 

Comparing ( |3.22| ) and ( |3.23j ) one concludes, that the derivatives in the 



r.h.s. of (|3.23|) act on the renormalized part of the W-loop as well as on 
the Z - factors. As a result a mixing between different terms appears and 
contour divergencies cannot be separated in a simple way in the case of FC, 
and the perturbative behaviour of FC is controlled by both ultraviolet and 
contour logarithmic divergencies. 



Direct perturbative analysis of the simplest nontrivial bilocal FC [|114|, 



115|| shows, that this is indeed the case. In a particular case where the 
reference point coincides with one of correlator arguments the perturbative 
expressions for the functions D and D\ at the next-to-leading order reads 



iV 2 -l 1 a s 
2 



N 2 - 1 1 a, r 1 
— L 



2 lt 2 Z A TT 



4 



(3.24) 



(3.25) 



A2JV C 4/ 3N C 24 3 

with L = log(e 27 /i 2 2; 2 /4), \x being a renormalization scale in the MS scheme 
and the first coefficient of the /9-function (3\ = (11N C — 2/)/6. 

Note, that at the tree level D-function receives zero contribution while 
£>i is: 

D?\z 2 ) = ^±4-,=---. (3-26) 
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Since the string tension is an integral over D(x) ( ft.24| ) tells that the 
0(g A ) perturbative contribution to D(x) produces nonzero string tension, 
which physically has no sense and should be cancelled by other terms of 
the cluster expansion. Such cancellation was indeed shown to occur |89| in 
perturbation theory, so string tension is exactly zero at any finite order of 
the perturbation theory in the gauge-invariant formalism considered, i.e. if 
the latter is developed in terms of FC. 



3.4 The string profile 

It was already mentioned above that the main phenomenon behind the 
physics of confinement is the formation of the string between the colored 
sources. It is this string which mediates attraction between quark and anti- 
quark with the force of roughly 15 tons. 

An immediate question arises: what does the string consist of? Or in 
other words; what is the distribution of fields inside the string? 

Before answering this question, two specifications are needed. First, 
speaking about field in the string, one should always compare it to the field 
in the vacuum, since only this difference makes the string. Second, in the 
QCD vacuum, which is stochastic by necessity as discussed above, the gauge 
invariant field distribution can be measured only with the help of field cor- 
relators. 

It is instructive to distinguish two general types of correlators in this 
respect: 

i) connected correlators - correlators of colored objects, which we dis- 
cussed above, where the trace operation appears only once (if at all) 

((Tr G^u, x )G p(7 (v, x ))) 

ii) disconnected correlators - correlators of two or more white objects, 

like 

«?:F%{x)TrF%{y))) 

In nonabelian theory, unlike in QED, those two correlators correspond to two 
very different probes, which can be made with fields in the string. In the 
first case one measures the color field distribution as it is, in the second one 
measures the "interaction" between two white objects: the string and the 
neutral plaquette, which can exchange only white objects - glueballs in case 
of gluodynamics. These approaches are complementary to each other and in 
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both cases one gets some gauge-invariant characteristics of the problem. In 
what follows we will compare the results of analysis performed with these 
two kinds of probes. 

We start with the case of the connected correlators. Using Monte-Carlo 



(MC) technique, in [g7|, H, H] the study was made of the spatial distribution 
of the components of the field strength tensor in presence of a static qq pair. 
Following |p8|| , we define: 

_ (Tr (W-$-P fw (S||,SL)$ + )) 

/v- {T[W) 1 l^-^j 

where PU is a W-loop, $ is a phase factor and P^ u is the plaquette, oriented 
in order to give the desired component of the field.[] The coordinates x\ux± 
measure the distance from the edge of the W-loop and from the plane defined 
by the loop, respectively. In the naive continuum limit a — > 0, 

Plxv ~ a 2 {F„ v ) qq (3.28) 

where (F q g) is the average value of the complicated Wilson contour with 
plaquette P^ v reduced to F^a 2 . In a 16 4 lattice was used, taking a 



8x8 W-loop and (3 = 2.50, which is inside the scaling window for the fields. 
Moving the plaquette in and outside the plane defined by the W-loop, one 
obtains a map of the spatial structure of the field as a function of x\\ and 
x±. Using a controlled cooling technique, one eliminates the short-range 



fluctuations |[23|| . The long-range non-perturbative effects survive longer to 
the cooling procedure, showing a plateau of 10-14 cooling steps, while the 
error becomes smaller. A similar behaviour has been observed for the string 
tension and the magnitude of string tension at the plateau coincides with 
the full uncooled value to a good approximation, which confirms the physical 



significance of string measurements, done in |27], |28], £9 . 

The cooling technique allows to disentangle the signal from the quantum 
noise with a relatively small statistics. The general patterns of the field 
configurations are briefly resumed in the following figures, while the complete 



results for the field-strength tensor are contained in Tables of ||2S| , |29j. Fig. 
5 represent a map of the spatial behaviour of the longitudinal component 
of the CE field. From Figs. 5, it is seen that the parallel electric field is 
squeezed in flux tube. 



1 Indices fi, v in (3.27) are not usual tensor indices but numerical arguments, informing 



about the chosen direction of the plaquette. 
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We denote direction along the qq axis x\\ = x±, while that of x± — x 2 , and 
the Euclidean temporal axis is x 4 . All the construction is taken at a fixed 
value of x 3 . 

Using the non-abelian Stokes theorem and the cluster expansion (see 
chapter 3) for one has 

Pixu(xi,x 2: Xi) = a 2 J da^x'^x'^k^ (3.29) 

where 

= ^Tr (E^, 0, x'^F^(x 1 ,x 2 , x 4 )$ + ) + (3.30) 

$ is the parallel transporter from the point (x[, 0, x' A ) to ( ) , and dots 

imply contribution of higher order cumulants containing additional powers 
of E x . 

We shall keep throughout this section only the lowest cumulants (con- 
taining lowest power of E\) and compare our prediction with the MC data 
of previous sections. The bilocal correlator A^ can be expressed in terms of 
two independent Lorentz scalar functions D{{x^ — x'^) 2 ), Di((x^ — x'^) 2 ) (see 
chapter 2) 

A u = D + D 1 + (hl + hl)^ (3.31) 
A 24 = (M 2 )^ , A M = (h 1 h 3 )^± (3.32) 

A 23 = 0, A 13 = (h 3 h 4 )^ ; A 12 = (h 2 h 4 )^± (3.33) 
Here h ll = {x - z%. 

Since all construction in Fig. 5 is at x 3 = x' 3 = we have h 3 = and hence 

P23 = P34 = Pl3 = (3.34) 

The only nonzero components are A 14 , A 24 and A 12 . 

When x 4 = ( and this is where measurements of p\ 2 have been done), 
P\ 2 vanishes because of antisymmetry of A 12 in h 4 . Hence only p 14 and p 24 
are nonzero, and only those have been measured to be nonzero. 

On the Fig. 5 (a,b,c) the field distributions (pu(xi, x 2 )) 2 are shown for 
different quark separations, R = T g ,5T g and 15T ff . The string profile, i.e. 
(Pi4(^ 2 )) 2 distribution in the middle of the string, for these quark separations 
is shown in Fig. 5 (d). One can see a clear string-like formation with the width 
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of 2.2T g . The mean value of the CE field in the middle of string does not 
depend on R for R greater then 5T g . 

To make comparison with data more quantitative, one can exploit the 
exponential form of D, D\ |24 



D^h 2 



£>i(0) expC-^l); D(h 2 ) 
1 



D(0) exp(-6\h\) 



(3.35) 



Dx(0) « -D(0); $i«<S 



Inserting this into Q3.29 ), we have 
^14(^1,^2; 0) = a 



rR 

2 / dx\ 
10 



T 
2 



dx a 



D{0) + D 1 {0)-^±^D 1 (0) 



-61, 



(3.36) 



with 



-x 



4 1 



/ll = X\ 



X 



1 J 



h 2 = h\ + h\ + x\; 



For P24 one obtains a similar expression, which vanishes for 22 = and in 



the middle of the string, which exactly corresponds to the data of [25 



Finally, we can make a detailed comparison of the prediction for p 14 in 
( |3.36| ) with present lattice data. One obtains a simple analytic result for 
p 14 (x2 = £j_) in case of a very long string. The transverse shape measured 
at the middle is given by 



Pl4 



2na 2 



D(0)(l + 6x 2 ) 



Di(0)-(5x 2 ) 2 



-SX2 



(3.37) 



P14 was calculated as a function of x±, x 2 from (|3.36|) keeping -Di(O) = ^D(0). 
One finds 5 _1 ~ 0.2 Fm , which is in good agreement with ||24|| . 
The results allow to predict all curves for other values of x» and x±: the 
agreement with the numerical results is very satisfactory p5| , |26 . 



The asymptotics of the string profile at large x± is shown to be exponen 
( P-37j ), just as the asymptotics of D,D ly measured in 



tial, see eq 

This in contrast to the behaviour inside the string, where the Gaussian-like 
flattening is observed. This effect is connected to the smearing effect due to 
integration in ( |3.36| ), yielding polynomial factors in ( |3.37|) . The radius of the 
string is close to 5 _1 = T g and hence gluon correlation length T g has another 
physical meaning - the thickness of the confining string. 



34 



One can define higher irreducible correlator 7° (the superscript c is writ- 
ten for connected, so the trace appears in ( |3.38|) only once) as follows 

T*a*{(FF') q ,-(FF') } (3.38) 

From ( |3.38| ) it is clear, that 7 C contains only double plaquette correlations. 
Most of the lattice data for 7 C are compatible with zero net effect, within 
two standard deviations (see [2~5| and recent analysis in [50]). 



The same line of reasoning can be elaborated in the case of qqq states, i.e. 
baryons. The reader is referred to the paper where the corresponding 



field profiles of the string were obtained. An interesting difference from the 
meson case considered so far is the vanishing of the average value of the 
CE field at the point of the string junction in the baryon, occuring due to 
symmetry reasons. This fact does not mean of course the vanishing of higher 
order correlators, measuring the fluctuations of the fields around the mean 
value (see [|116|| ). 

For the double correlator *j c one can use the non-abelian Stokes theo- 
rem and the cluster expansion, to connect 7 C to the quartic correlator G 4 
and hence to estimate it from the lattice data. Representing the coordinate 
dependence of G4 as an exponent similarly to G2 = D(h 2 ) with the same 5, 
one gets the dimensionless ratio as 

G 4 (0) 7 C (0) ^ 2 _ 3 



(g 2 (o)) 2 (p(o)y 

From the point of view of the cluster expansion convergence, the expansion 
parameter is of the order G2(0)T^ <C 1, and one finds that G 4 (0)Tg ~ (2 — 3) ■ 
(G2(0)Tg) 2 in agreement with the idea of the bilocal correlator dominance. 

We turn now to the case of the disconnected correlators. The advan- 
tage of this probe is that it measures the field configurations induced by 
a static quark-ant iquark pair which are gauge invariant by itself, i.e. in 
the simplest case the traces of the squares of the CE and CM field compo- 
nents: J2a {Ei) 2 , J2a {Bi) 2 - These expressions are directly related to phys- 
ical quantities as for instance the energy density. We are interested in the 
difference between the squared field strength of the situation with a static 
quark- ant iquark pair and the vacuum expectation value without such pair. 
This difference which we denote by AF 2 ^ (x) is given by: 



Ap2 f v _ 4 (Wg] TrP aP (x)) - (TrW[C])(TrP aP (x)) 
M a p{x) - a4 (Tr W[C}) ' ( °-°" ,) 
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Here Tr P a p is the W-loop of a small square in the a — (3 plane. From it 
we obtain the desired gauge invariant squared field strengths: 



Tr P, 



a/3 



X 



(3.40) 



We see that we have to deal here with a more complex problem than 
in the previous case, namely the evaluation of the expectation value of two 
loops. Therefore the factorization hypothesis of the Gaussian model is tested 
here much more severely: vacuum expectation values contain at least four 
gluon field tensors which have to be reduced to two-point correlators by the 
factorization assumption discussed above. 

We again apply the non-abelian Stokes theorem to the line integrals in the 
loops and transform then in surface integrals as discussed before. Expanding 
the exponentials in eq. (|3.39|) we obtain: 



x 



4 1 



surfaceordered 



£< 



, n=l 



da 



w 



da™ v Tr \t ai 



dolvdo p p(J x 



-i 



-Tr 



t a t b 



V p\V\ ' ' ' 1 /J, n u n 



surface— ordered 



£(- 

n=l 



da w 



■ ■ Tr [t° 



pa pb \ 
fiv pal 

••f" 



■t an ] X 



d°l v do p pa x 



2! 



-Tr 



I p a i . . . p a " \ I T? a p b \ 

X- 1 pii/i ^vMnl \ pv pal 



(3.41) 



The surfaces Sw and Sp with the elements da w and da p are bordered by 
the loop of the static quarks and the plaquette respectively and contain both 
the reference point y, the choice of these surfaces will be discussed later. The 
indices of the surface elements are restricted to fi < v . 

Using the factorization prescription given in the previous section in and 
subtracting out the terms where the fields on Sw and Sp factorize separately 
we obtain: 



Ai^ (x) 



4 1 



a 4 (W) ( 
Tr [t ai ■■■t a2n 



surfaceordered 



\2n 



n=l 



dcr™,., ■ ■ ■ d(T^ n ,.„ X 

P\V\ P2n,V2n 



(3.42) 



\^ p ai v ai ' 

pairs 
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Of the 2n field strengths running over the surface Sw two are correlated 
with a F^ u on the pyramid of the plaquette and the other (2n — 2) among 
one another. We now use the prescription developed in the previous section 
and neglect (in any order of the expansion of the exponentials) all terms 
which are not fully ordered, i.e. which are suppressed by at least a factor 
Tg/L. For all fully ordered expressions the traces in ( |3.42j ) are identical and 



one obtains finally after a lengthy calculation 45 : 



1 



a 4 2iV c (iV c 2 - 1) 



1 2 



(3.43) 



We note that in the integral one surface element is connected to the 
surface Sw (indicated by the index W) and one to the plaquette oriented in 
a, /3-direction (indicated by P) As noted above the choice of the surfaces is 
somewhat arbitrary and depends on the way by which we transport colour 
flux from the point x to the reference point y. Choosing a straight line it 
is reasonable to put the reference point between the plaquette and the W- 
loop. In this case the surfaces are the sloping sides of pyramides with the 
reference point as apex. It was checked in [117|, [45 , that different choices 



of the surfaces leads only to moderate changes in the results. A reasonable 
strategy seems to be to choose the reference point in such a way that the 
total resulting surface is minimal. 

For the evaluation of ( |3.43| ) we refer to the paper [|117|| , p5| and quote only 
the main results in below. First of all, it can be seen by symmetry arguments 
that for that case the difference of the square of the magnetic field strengths, 
i.e. the plaquette P a p{x) with no time (24) component, vanishes identically. 
This means that the CM background field is not affected by the static color 
charges. 

The electric field perpendicular to the W-loop is also practically not 
affected but only the difference of the squared electric field parallel to the 
loop is. In both cases the squared electric field difference is negative, i.e. the 
presence of the static source diminishes the vacuum fluctuations 

In Figs. 6, 7 we display the value of the squared field strength parallel to the 



spatial loop extension AF ; 
coordinates r = 



(x) = -(E^x 3 ,r)) 



qq— vacuum 



as a function of the 



x\ + x\ and X3 (the results are rotational invariant around 
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the a^-axes) for different spatial separations Rw (notice, that a = 0.35 fm 
on these pictures is a chosen scale of distance analogous to the correlation 
length T g and has nothing to do with the lattice spacing a used above). 

The squared field strength reaches its saturation value AF 3 2 4 (x) which 
approximately equals to 14MeV/Fm 3 for a spatial extension of the W-loop 
Rw ~ 4T 9 . The transversal extension, defined by 



Jdrrr 2 AF| 4 (0, r) 
\ JdrrAFi 4 {0,r) ' 1 } 

is practically independent of Rw and about 1.8 times the correlation length 
T g . As mentioned above, the quantity AF 3 2 4 (x) is the difference between 
the squared field strength in the presence of a static gg-pair and the non- 
perturbative vacuum field strength. One can construct several quantities of 



interest from it. For example, the figure 10 from the reference ||117|| displays 
the unsubtracted squared field tensor of the color source in units of the vac- 
uum gluon condensate (F 2 ), while the figures 11-13 from the same reference 
shows the squared field strength including the perturbative, i.e. Coulomb 
contribution. 

The result is therefore very satisfying for a physical intuition: a CE flux 
tube is formed between the quark-antiquark pair whose energy content in- 
creases linearly with the separation of the pair. One may even compare the 
energy content of the flux tube formed in that way with the quark antiquark 
potential calculated with the help of a single W-loop. In order to do that we 
discuss two low energy theorems in the next subsection. 

The analysis of higher cumulants contributions for the disconnected cor- 



relators was performed in ||118j| . It was shown, how quartic cumulant can 
modify the profile of the string and, in particular, lead to (small) nonzero 
components of the magnetic fields averages. 

The pictures displayed on the Fig.5 and Figs. 6, 7 are rather similar, despite 
they were obtained by use of different probes, and both are in good agreement 
with the lattice results for the string profile. We found it instructive to 
demonstrate how FC formalism works in these two different settings: in the 
first case the mean value of the CE field in the presence of the static pair 
was calculated in the Gaussian approximation and its squared average was 
shown on the Fig.5; in the second case the correlation of two gauge-invariant 
object was taken and reduced in the Gaussian approximation to the integral 
of the square of the two-point correlator again. From general point of view 
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the typical correlation length, characterizing the decrease of the fields in the 
case of connected correlators is related to the mass of the gluelumps (and to 
T g x in the Gaussian approximation), while that of the second case represent 
somehow the glueball mass, since it corresponds to the interaction of two 
white objects. 

Summarizing, all the results discussed so far, together with the arguments 
on Casimir scaling supports the fundament ality of the bilocal correlator, 
which defines the nonperturbative dynamics of confinement. On the contrary, 
there are no data known to the authors up to now, which would indicate a 
significant role played by the highest field strength correlators. 

3.5 Low energy theorems and consistency of the model 

In this subsection we consider two low energy theorems and apply them to 
the results of the previous sections in order to check the consistency of the 
model. We follow here closely the line of reasoning given in [117]. Let us 
consider the vacuum expectation value of a W-loop: 

(W[C\) = ^(TrPexpiiJ Apdxp]) (3.45) 

Here C denotes a closed loop. Differentiating log(W[C]) with respect to — ^ 
one obtains: 

Tr Pexp[i J A^dx^S - (J d 4 xTr [F„ v (x)F^(x)]j (3.46) 

If the loop C is a rectangle with spatial extension —R/2 < x% < R/2 and 
temporal extension —T/2 < £4 < T/2 the right-hand side of ( |3.46| ) can be 
simplified further in the limit of large T. Let the length T g be of the order of 
the correlation length of the gluon field strengths. Then for \x±\ > T/2 + T g 
the r.h.s. of ( |3.46|) is zero and for \x±\ < T/2 — T g independent of x 4 . This 



gives in the limit of large T: 

8 ^l0g[(W/[C])] = / I [p {x) p {x)] \ . (3 47) 

oa s0 \J I r 
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Here and in the following x 4 is zero and the expectation value {..)r shall 
denote the expectation value in the presence of a static quark-antiquark pair 
at distance R, with the static sources in the fundamental representation, and 
the expectation value in absence of the sources has been subtracted. 

The potential V(R) is defined by (|3.15|) and thus one obtains the low 
energy theorem: 

~ 4 ™*°?a~7 = \ (I dV& FrW^M) B ( 3 - 48 ) 

In ( |3.48|) self energies are understood to be subtracted. 

Up to now we have used the bare coupling g and bare fields. If we use 



the background gauge fixing [121, 122, 123 1 we have 



g = Z g l g Q , A, = Z-^Af , Z g Z 1 / 2 = 1 , A, = g Af = gA, (3.49) 

where g and A^ are renormalized quantities. The renormalized squared field 
strength tensor is given by: 

Tr [F, v {x)F, v {x)\ = Z^Tr [F$(x)F<${x)] (3.50) 



where Z F i was calculated [|124j| in background gauge with Landau gauge 
fixing to be Z F 2 = 1 + \go-^ log Z and thus we obtain 

|^ = ZZ F 2 . (3.51) 
da s0 

We use ( p.49| )-( p.51| ) to express ( ^.48| ) in terms of renormalized quantities: 
dV(R) 1 



a* 



da. 



J d 3 xTr [F^{x)F^{x)}^ (3.52) 



In order to evaluate the left-hand side of ( |3.52| ) we use standard renormal- 
ization group arguments: the only scales are the renormalization scale \x and 
R, thus the potential V = V(R, /i, as) satisfies on dimensional grounds: 

{ R M ~ V{R ' /i ' as) = ~ V{R ' ^ as) (3 ' 53) 

But for V as a physical quantity the total derivative with respect to fi must 
vanish: 
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d ( d d \ 

li-^-V(R, n, a s ) = I /I— + P( a s)-Q^J v ( R i ^> a *) = 

where (3(a s ) = n±a s {fj,). We get from (|33^ ) and ( ggg) : 
dV(R) _ 1( dV(R)\ 

and finally one obtains with (|3.52D the low energy theorem 



V(R) + R^^) =l£(J d " xTT \F^)F^)\)r (3-56) 



or equivalently 



(3.54) 



(3.55) 



(Jd 3 x(E(x) 2 + B(x) 2 )) R (3.57) 



The ( p.57| ) is for the case of a linear potential just one of the many low energy 



theorems derived in [35]. But there renormalization was only discussed to 
leading order. 

The relation between the potential V(R) and the energy density T 00 (x) 

is 



V(R) = / d 3 x(T 00 (x)) J 



^(Jd'x (-E( 



+ B(xr)) R 



(3.58) 



Note that we are in a Euclidean field theory, hence the minus sign of the 
electric field. Moreover on the r.h.s. of (|3.57|) and (|3.58|) we have renormalized 
composite operators. 

In the model under consideration where the CM and the transversal com- 
ponent of the electric field vanishes the potential must be proportional to the 
spatial integral over E 2 . In view of the approximations made concerning path 
ordering as well as the arbitrarines of the choice of the reference point this 
is a very valuable consistency test (see details in ||117|| ). 

In lattice QCD (|3.57| ) and (|3.58|) are known as Michael's sum-rules ||125 



But in the derivation of the action sum-rule ( |3.57| ) scaling of V(R) with R 
was not taken properly into account and hence the second term on the l.h.s. 
of ( p.57|) is missing in ||125 . 
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First consider the region where the potential is linear. Then ( p.57| ) reads: 



2 a, 



■(|d 3 a; (E(x) 2 + J B(x) 2 )) /? = 2\/(^) 



Since in the model of the stochastic vacuum the static sources do not modify 
the CM field, i.e. (B 2 )r = , we have consistency between ( |3.58| ) and (|3.59|) 
only if 

(3/a s = -2 (3.60) 

From the energy sum-rule ( |3.58| ) we have already obtained a s = 0.57, thus 
consistency of the MSV requires that (|3.60|) should be satisfied for this value 
of a s . 

The /5-function can be expanded in a power series in a s : 



(5(a s 



A) 2 A 3 4 . 

a or or + 

2tt s 4tt 2 s 64tt 3 s 



(3.61) 



The n-loop expansion gives the series up to a° +1 . The first two terms 
of the power series are scheme independent, the numerical value of the term 
proportional to a A s has been evaluated ||126|| in MS. For pure gauge SU(3) 
theory we have M: 



/3o = 11, Pi = 51, /3 2I7 -= 2857 



(3.62) 



We see that the condition ( |3.60| ) yields a s = 1.14 using (3 on the one 
loop level, a s = 0.74 on the two loop level and a s = 0.64 on the three loop 
level in MS. This value is already very close to the value of a s = 0.57 
determined from the flux tube calculation using (|3.58|) . If we use for f3 the 
series (|3.61| ) truncated at order and impose condition Q3.60D for a s = 0.57 
and determine the coefficient 0% accordingly we find @2 = 6250 i.e. about 
twice the MS value. 

It is also easy to see from the above that for distances R/T g < 1, where 
T g fa 0.3 fm the MSV does no longer yield a linearly rising potential and 
hence we have 



V(R) + R 



dV{R) 
OR 



+ 2V(R) 



R/T g < 1 



Therefore the conditions Q3.58] ) and ( p.59| ) cannot be fulfilled simultaneously 
for (B 2 )n = 0. This does by no means speak against the model, since in the 
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derivation of the results it was stressed that the spatial extension of the loop 
had to be large as compared to the correlation length of the field strengths 
in order to justify the evaluation of E 2 and B 2 . 

Notice, that no use of the equations of motion has been made explicitly. 
The factorization hypothesis (Gaussian measure) enters crucially in the eval- 
uation of the squared color fields. In view of the drastic approximations made 
it is gratifying that relation ( |3.60|) , which is a highly non-trivial consequence 



of the dynamics of the system, is so well fulfilled. 

The scale /i at which the MSV works was found to be the one where 
a s (/z) = 0.57 by using the energy sum-rule. With the help of the action sum- 
rule we now checked the consistency of the MSV at this scale and obtained 
an effective /^-function. 

As can be seen from (|3.57|) and (|3.58|) the squared field strengths B 2 and 



E 2 depend on the renormalization scale fi and hence on a s (/x). We may use 
relations (|3.58| ) and ( |3.59| ) in order to predict the ratio of 



fd*x(B(x) 2 ) R _2 + (3/a s 



Jd 3 x(E(x) 2 ) R 2-p/a, 

as a function of the strong coupling a s . This result can be checked in lattice 
gauge calculations. 
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4 Heavy quark dynamics and field correlators 



In this chapter we apply the MFC to heavy quarkonia having in mind first 
of all bottomonium and charmonium and, qualitatively, ss states. For these 
systems one should take into account both perturbative and nonperturbative 
interactions, since the size of ground states is around 0.2 and 0.4 Fm for 
bottomonium and charmonium respectively, and the color Coulomb interac- 
tion is dominant below 0.3 Fm. As we shall see in this chapter the MFC 
in general and the Gaussian stochastic approximation in particular are very 
well adapted to the description of heavy quarkonia. 

Actually MFC is here the only analytic method which takes into account 
both perturbative and nonperturbative interaction at all distances. The first 
account for the effects of finite correlation length T g was made in Wl 



in 



attempt to correct the unphysical n 6 -growth of NP energy shift obtained in 
]60| , [61]. For states of large temporal extent (here belong all states of cc 



and ss) the potential picture emerges from MFC with potentials expressed 
through the only two correlators D and D\. For states of small size the NP 



configurations can be accounted for in form of condensates, (see Voloshin [60 



and Leutwyler [pl||). The general situation was discussed in |35], |39 



Various comparisons and cross-checks can be made with the MFC results. 
First, one compares with numerous experimental data and finds in most cases 
good agreement; second, lattice MC data for scalar and spin-dependent po- 
tentials are available; third, phenomenological potential models in conjunc- 
tion with the 0{a 2 s ) perturbative results are checked versus MFC predictions. 
The reader is also referred to the review ||127|| where different aspects of the 
heavy quark dynamics in the confining vacuum are extensively discussed. 

The material of present chapter is based mostly on papers |35|, [SB] and 
lectures |59| and the paper ||128| . 

The structure of the chapter is as follows. In section 4.1 the general 
formalism of the Feynman-Schwinger representation (FSR) for the quark- 
antiquark Green's function is presented and the nonrelativistic approxima- 
tion is deduced. In section 4.2 the NP and perturbative parts of static poten- 
tial are derived from the MVC. In section 4.3 the spin-dependent potentials 
are obtained from MVC and compared to the lattice data. Finally, in section 
4.4 results of numerical calculations and general discussion are given. 
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4.1 The general formalism for the quark— ant iquark Green's 
function 

We start with the quark Green's function in the form of the proper-time 
and path integral (the Feynman-Schwinger representation [35, 36. 37, 38, 

ds / Dze~ K T(x,y) (4.1) 



(4.2) 



S(x, y) = i(m — D) 
where T contains spin insertions into parallel transporter 

rx rs 

Y(x,y) = P A Pf exp[z / A^dz fl + dr gEF(z(r))), 

Jy JO 



K = m s + 



A Jo 



z^dr, 



and double ordering in and F^ u is implied by operators Pa, Pf- 
We have also introduced the 4x4 matrix in Dirac indices 



P>i Ej 



(4.3) 



Neglecting spins one has instead of ( 4.2 



rx 

T(x,y) — ► $(x,y) = Pexp i / A^dz^ 

Jy 



(4.4) 



In terms of (qq) Green's functions ([4.1|) and initial and final state matrices 
Ti, Tf ( such that q(x)T /$(x, x)q{x) is the final qq state) the total relativistic 
gauge-invariant qq Green's function in the quenched approximation is 



G(xx,yy) = (Tr (T f Si(x, y)T t ^(y, y)S 2 (y, x)$(x,x))) 

- (Tr (TfS^x, x)$(x, x))Tr (T t S 2 (y, y)<$>(y, y))) 
The angular brackets in ( |4.5| ) imply averaging over gluon field A^. 



(4.5) 



The second term on the r.h.s. of (O) represents a mixing of a flavour 
singlet quarkonium with two- or three-gluon glueball. Except for an acci- 
dental degeneration of masses, this term is a small correction and will be 
omitted in what follows. 

Since we are interested primarily in heavy quarkonia, it is reasonable to 
do a systematic nonrelativistic approximation. To this end we introduce the 
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real evolution parameter (time) t instead of the proper time r in K, (r in K) 
and dynamical mass parameters /i, p 

dt dt r s 9/ , , r T , ( dz u (t)\ 2 

Here we have denoted T = (X4 + X4) /2. Nonrelativistic approximation means, 
that one writes for 24 (t), 24 (t) 

Zi (t) = t + ((t), 2/i 1 = -; (4.7) 

si 

z 4 (t)=t + ((t), 2/i 2 = - 

and expands in fluctuations £, £, which are 0(-j=) . Note that the integration 
in ds\ds2 goes over into d\i\d\i 2 . Physically expansion ( |4.7| ) means that we 
neglect in lowest approximation trajectories with backtracking of 24, 24, i.e. 
neglect qq pair creation. One can persuade oneself that insertion of ([4-7]) into 
i^, .K" allows to determine fii,fi 2 horn the extremum in iT, ^ to be 

fix = m 1 + 0(l/mi), ^ 2 = m 2 + 0(l/m 2 ) (4.8) 

and one can further make a systematic expansion in powers 1/mj. (For details 
see [j35]). At least to lowest orders in 1/m* this procedure is equivalent to 
the standard (gauge-noninvariant) nonrelativistic expansion. 
Let us keep the leading term of this expansion 

G(xx, yy) = 4m 1 m 2 e- {mi+m2)T J D 3 zD 3 ze~ Kl - K2 (W(C)) (4.9) 

where K x = ?f J T ' zf(t)dt, K 2 = ^ ^{t)dt and (W(C)) is the W-loop 
operator with closed contour C comprising q and q paths, and initial and 
final state parallel transporters $(x,x) and $(y,y). 

The representation fl4.9Q is our main object of study in next section. 



4.2 Perturbative and nonperturbative static potentials 

The technical tool we are using is the background perturbation theory. We 
represent the total gluon field as 

A lt = B^ + a lt (4.10) 
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where _B M is the NP background, while is perturbative fluctuation. The 
principle of separation in ( |4.10|) is immaterial for our purposes (for more 
discussion of this point see [ |129|| ). 



For static charges trajectories are fixed straight lines and (W(C)) in ( f4.9|) 
factorizes out defining the potential between static charges V(R) 



(W(C)) = exp(-V(R)T) 



(4.11) 



where R is the distance between charges and T is the time distance in the 
rectangular loop; C = R x T. 

In the approximation, when only bilocal correlators are kept, V(R) is 
given as |08l |39| 



V (R) = C D 



a s {R) 
R 



rR roo 

2R dX duD(X,u)+ 
Jo Jo 



+ / XdX / dv(-2D(\,v) + D 1 (\,v)) 
Jo Jo 



(4.12) 



In addition one also obtains radiative corrections due to transverse gluon 
exchange, which have been computed earlier ||130 . 
Thus one has the resulting potential 



V{R) = V (R) + AV(R) 



(4.13) 



where Vq(R) is defined in ( 4.12 ), and AV(R) is the sum of higher-order 
terms, namely 

AV(R) = AV^R) + AV 2 (R) + AV 3 (R), 



where AVi(R) refers to 0{a 2 s ) perturbative contributions ||131| , |132|| , AV2(R) 
is interference term [ |133|| which is 0(a s ) and AVs(R) refers to contribution 
of higher FC. 



4.3 Spin— dependent potentials 

As was shown in the previous section, the leading term in the heavy mass 
expansion is the nonrelativistic form containing the W-loop with the contour 
C, formed by straight-line paths without spin insertions (S-F), Eq. ( |4.9| ). To 



keep spin dependent terms we shall use (4.2), (4.5) and derive from it the 
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nonrelativistic expression in the limit when m; — > oo and also \x — y\ <C 
\xq ~ 2/o I - 

Now one can use the standard definitions of Eichten, Feinberg and Gromes 
(EFG) |3|, [135| for the SD potentials 

(JiLx (X2L2 



Vsd(R) 



Am 2 



1_ de_ + 2dV 1 (R) 
R dR RdR 



+ 



+ 



v Am{ <±,,o 2 

azLt-aiLi l dV 2 {R) | a^V^R) | (SoiRagR^^gg^g) (4 M) 
2m\m 2 R dR Ylm\m 2 V2m\m 2 R 2 
and identify Vi(.R) by expanding the exponents of (Si-F) in (|4.2| ) to the first 
order and keeping m~ 2 terms as compared to the leading spin-independent 
(static) contribution. 

In Gaussian (bilocal) approximation one obtains (we denote quantities 
taken in Gaussian approximation by primes) 



1 dV 2 ' 
R dR 



1 dVi 
R~dR 



dv 



dv 



RdX 

"r 



R. 



R \d\ 



R 2 



D(X, v 
D(X, v) + D 1 (X, v) + \ 



2 dDi 
d\ 2 



1 d,p' 1 r x / // r" 
-RlR = -RL d »{l D ^ ) + L 



R 



0R 2 



Vi = r dv \3D{R, v) + SD^R, v) + 2R 
J —00 V 



,dD x 



(4.15) 

(4.16) 
(4.17) 

(4.18) 
(4.19) 



dR 2 J 

It is very simple to check, as it was done in [^] , that the Gromes relation 

V{ - V 2 ' + e' = (4.20) 

is identically satisfied in bilocal approximations for e', V(, V 2 given by ( |4.15| ), 
( 1.16| ) and ( 4.17 ) for the arbitrary choice of D and Di falling off at large 
distances. 

One can also check that (|4.17|) is the derivative of the NP part of the 
scalar potential fl4.12| ) (without Coulomb part). 

It is interesting that fl4.15|) and ( |4.16| ) yield correct results in lowest order 
of perturbation theory, namely V[ has no contribution and is equal to 
when the perturbative value of Di, D'{ ert {x) = 16a s /(37ra; 4 ), is used. 
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4.4 Charmonium and bottomonium spectra. 



In previous sections the scalar and spin-dependent potentials have been ex- 
pressed through the Gaussian FC, functions D and D\. The use of these 
potentials in the Schrodinger equation allows to calculate the spectra of char- 
monium and bottomonium and compare results with experimental data. The 
contents of this section is mostly based on |136|. 

The static potential according to can be written as 



V (r) = V per \r) + e{r) 



(4.21) 



The NP component of D, D\ are taken from lattice measurements [g4J] in the 
form (see Section 2) 



D{x) = D(0)e- {xl/Ts , Di(x) = £>i(0)e~ 



\*\/T g 



(4.22) 



where D(0) « 3£>i(0). 

Parameters used above for computations and corresponding to lattice 
results |22|, 1241 are 



a = 0.185 GeV 2 , a x 



-D 1 (0)T ( 

71 



t: 



1 GeV, 



4.85 GeV, 



0.0064 GeV 2 , GeV 
1.45 GeV 



m 



(4.23) 

The perturbative part of potential takes into account the asymptotic free- 
dom and freezing at large distances with a s = a™ ax = 0.5 32, 33]. 

Results of computations for the parameters of Eq. ( |4.23| ) are listed below 
in the Tables 1-2 together with experimental data and calculations with the 
Cornell potential | |137| . 
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Table 1 



Masses of bb for FC (4.22) in comparison with experiment and those of 



Cornell potential [137 



State 


Exper. 


FC 


Cornell 


7(15) 


9460 


9460 


" fit" 


7(25) 


10023 


10023 


10052 


7(35) 


10355 


10347 


10400 


Xb(lP) 


9900 


9918 


9960 


Xb(2P) 


10260 


10250 


10314 


7(25) -7(15) 


563 


563 


592 


7(35) -7(25) 


332 


324 


345 


7 (2S) - X6 (1P) 


123 


105 


93 


Xb (2P)- X b(lP) 


359 


332 


354 



Table 2 



Masses of cc for FC ([4.22 ) v s experiment 



State 


Exper. 


FC 


Cornell 




3097 


3096 


3097 


*(2S) 


3686 


3682 


3688 


*(3S) 


4040 


4128 


4111 


Xco 9 (lP) 


3525 


3516 


3525 


*(1D) 


3770 


3800 


3810 


^(25) - J /if;' 


589 


586 


591 


^{2S)-Xco 9 {lP) 


161 


166 


163 


V(1D) - J If 


673 


705 


713 



Results listed in Tables 1-2 demonstrate several successes of MFC in com- 
parison with the standard potential models like Cornell. Namely, i) it is 
possible to describe bottomonium spectrum with the current 6-quark mass, 
mj = 4.8 GeV, while in Cornell potential an effective mass rrib = 5.17 GeV 
is used; ii) a good description of leptonic width is obtained (not shown in 
Tables 1,2, see ||136| ) ; iii) the overal agreement of levels with experiment is 
rather good and certainly better than for the Cornell potential. 
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The fine structure of the spectrum is much more delicate issue. In particu- 
lar, here one should use the 0(a 2 s ) contributions to LS and scalar interaction. 
The full analysis done in ||138|| demonstrates a good description of charmo- 
nium and also checks the consistency of the model with the lattice data on 
spin-dependent potentials ||139| . 
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5 Field correlators at T > and deconfine- 
ment transition 



The finite temperature QCD provides a unique insight into the structure of 
the QCD vacuum, in particular it gives an important information about the 
mechanism of confinement. It is also of more practical interest, since the 
nature of the deconfinement transition has its bearing on the cosmology and 
the hot QCD plasma can possibly be tested in heavy ion collisions. 

From theoretical point of view the hot QCD is a unique theoretical labo- 
ratory where both perturbative (P) and nonperturbative (NP) methods can 
be applied in different temperature regimes. 

It is believed that the perturbative QCD is applicable in the deconfined 
phase at large enough temperatures T, where the effective coupling constant 
g(T) is small [ 140 , 141 ], while at small T (in the confined phase) the NP 
effects instead are most important. However even at large T the physics 
is not that simple: some effects, like screening (electric gluon mass), need a 
resummation of the perturbative series ||142| , |143|| , while the effects connected 
with the magnetic gluon mass demonstrate the infrared divergence of the 
For the role of hadron degrees of freedom in the transition see 



series 


144 


e.g. 


145 





During last years there appeared a lot of lattice data which point at the 
NP character of dynamics above T c . Here belong: i) area law of spatial 
W-loops ||146 , P-4711 , ii) screening masses of mesons and baryons [|148| , |i~4l 



and glueballs ||150|| , iii) temperature dependence of Polyakov-line correlators 



151 



In addition, behaviour of e — 3p above T c has a bump incompatible with 
the simple quark-gluon gas picture [ |152j| . 

Thus the inclusion of NP configurations into analysis at T > and also at 
T > T c is necessary. This was done using background field method in ||153 



Besides defining NP dynamics, the NP configurations also ensure freezing of 
a s and stabilize perturbative series |153|. To describe the phase transition a 



simple choice of deconfined phase was suggested ||154|| where all NP CM con- 
figurations are kept intact as in the confined phase, whereas CE correlators 
responsible for confinement, vanish. 

Recently this picture of phase transition has been tested and confirmed 
in the lattice calculations [25, ^(J of electric and magnetic correlators. It was 
shown in |27|, |23|, |2P| that the confining part of the electric correlator D E 
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decreases steeply right above T c , while the deconfining electric part, Df , and 
magnetic part D B , Df, stay approximately constant, see Figs. 2, 3. Thus one 
has a consistent mechanism describing the structure of the QCD vacuum at 
T > and the phase transition at T = T c . This picture together with the 
background perturbation theory forms a basis of quantitative calculations, 
where field correlators (condensates) are used as the NP input [154 . 



The plan of the Section is as follows. In subsection 5.1 the modified back- 



ground field formalism ||153|| is presented, based on the familiar background 
field method ||121 , |122 , |123|| for T > and incorporating the 't Hooft's iden- 
tity for integration over quantum fluctuations as well as over the background 
fields. The temperature phase transition is discussed in subsection 5.2 and 
resulting predictions for T c are compared with lattice data. In subsection 5.3 
the spacial W-loops are computed both below and above T c . 

5.1 Background field formalism 

We derive here basic formulas for the partition function, free energy and 



Green's function in the NP background formalism at T > ||154|| . The total 



gluon field is split into a perturbative part and NP background 

A„ = B„ + a„ (5.1) 

where both and subject to periodic boundary conditions. The 

principle of this separation is immaterial for our purposes here, and one can 
average over fields B^ and a M independently 

DA, exp(-S W ) = SDB^fDa^M-SiB + aH ^ 

= ((exp(-S(B + o))}„} B 

with the weight i](B). In our case we choose t](B) to fix field correlators and 
string tension at their observed values. 

The partition function can be written as 

Z(y,T,ii = 0) = {Z(B)) B 

where 

Z(B) = N J Da^ exp(-^ dr J d 3 x L(x,r)) (5.3) 
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where denotes all set of fields a M , N is a normalization constant, and 

the brackets (■■■) B means some averaging over (nonperturbative) background 
fields with the weight f]{B), as in ( |5.2| ). Explicitly 



8 ® 

5(5 + a) = [ dr [ d 3 xY I 



where 



In = \(F; v (B)) 2 ; I 2 = \alW; b u al 
L 3 = & a (D 2 (B)) ab e b ; U = -ig® a {D^ a^ ah e h (5.4) 

Ik = — ( 

5 2 

I 7 = -a v D^B)F^{B)- I 8 = *+(m + D(B + a))tf 



-(^(^(S)^) 2 ; I 6 = I int (a 3 , a 4 ) 



Notice, that we have redefined the fields according to — > gA^ in this 
section. 

In ( |5~5] ) 0, 9 are ghost fields, a is gauge-fixing constant, L 6 contains 
3- and 4-gluon vertices. For the fields B^ satisfying classical equations of 
motion Ij = 0. 

The inverse gluon propagator in the background gauge is 

W£ = -D\B) ab ■ 6, u - 2gF; u {B)r b (5.5) 

where 

(D x ) ca = d x 5 ca - igT b ca B b x = d x S ca - gf bca B b x (5.6) 

We consider first the case of pure gluodinamics, L 8 = 0. Integration over 
ghost and gluon degrees of freedom in (|5.3|) yields 

Z(S) = N'{VetW{B)Xet [Bet(-D,(B)D,(B + a))] a= ^_ x 

x{l + £%fa = ^)expf-V-V), (5.7) 



J„= D^B)F^{B) 



One can consider strong background fields, so that g-B^ is large (as com- 
pared to Aq CD ), while a s = |- in that strong background is small at all 
distances |JT], |32|. In this case ( |5.7[ ) is a perturbative sum in powers of <7 n , 
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arising from expansion in (ga^) n . In what follows we shall discuss the Feyn- 
man graphs for the free energy F(T), connected to Z(B) via 

F{T) = -T \og(Z(B)) B (5.8) 

As will be seen, the lowest order graphs already contain a nontrivial 
dynamical mechanism for the deconfmement transition. To the lowest order 
in ga^ the partition function and free energy are 



Z = (exp(-F Q (£)/T)) 



B 



F (B)/T = - log DetW - log Det(-D 2 (B)) = 

= Tr fat)j(-\e- tw + e tD2 ^) (5.9) 

where W = —D 2 {B) — 2gF and D 2 {B) is the inverse gluon and ghost prop- 
agator respectively, ((t) is a regularizing factor \\L54\\ . 

The ghost propagator can be written as [ |154j , |155| , |156|| 

(-D 2 )-^ = (x\ j o ds exp(sD 2 (B))\y) = J Q ds(Dz)™ y e- K <5>(x,y) (5.10) 

where winding path integral is introduced | 154 | 



(Dz) w = lim TT - C(m) T I dp e HY,«™)-(*-y)-^) ( 5 n) 



m=1 (4vre)2 n= ^ ±l J (27T 



with (3 = 1/T. For the gluon propagator an analogous expression holds true, 
expect that in (|5.5|) one should insert gluon spin factor P^exp 2gF inside 
<&(x,y). For a quark propagator the sum over windings acquires the factor 
(— l) n and quark spin factor is exp (gcr^F^) ||154 . 



5.2 Temperature phase transition in QCD 

We are now in position to make expansion of Z and F in powers of ga^ (i.e. 
perturbative expansion in a s ), and the leading-nonperturbative term Zq,Fq 
- can be represented as a sum of contributions with different powers of iV c 
and we systematically will keep the leading terms 0(N 2 ), 0(N C ) and 0(N®). 

To describe the temperature phase transition one should specify phases 
and compute free energy. For the confining phase to lowest order in a s 
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free energy is given by ( |5.3| ) plus contribution of energy density e at zero 
temperature 

F{1) = sV 3 - ^ 3 T 4 - T£ V ^f 2 ? ' 2 ex P (-^) + 0{1/N C ) (5.12) 

where £ is defined by scale anomaly 

1 1 n 

e ^ - l -N c ^-{{F; v {B)Y) (5.13) 

and the next terms in ( |5.12|) correspond to the contribution of mesons (we 
keep only pion gas) and glueballs. Note that e = 0(N%) while two other 
terms in ( pl2|) are 0(N®). 

For the second phase (to be the high temperature phase) we had made 
an assumption (which was proved on the lattice JST], |2B|| ) that all CM field 
correlators are the same as in the first phase, while all CE fields vanish. Since 
at T = CM correlators and CE correlators are equal due to the Euclidean 
0(4) invariance, one has 

((F; U (B)) 2 ) = {{F^f) el + ((Fp*) magn , (F 2 ) magn = (F'U (5.14) 

The string tension a which characterizes confinement is due to the electric 
fields |L6|, [L7|], e.g. in the plane (i4) one has 

o = (J E = g — J J d 2 x(Ti Ei(x)$(x, 0)^(0)$(0, x)) + ... (5.15) 

where dots imply higher order terms in (notice, that we have redefined 
the fields A — > gA in this chapter with respect to the notations adopted in 
Section 3.) 

Vanishing of oe liberates gluons and quarks, which will contribute to 
the free energy in the deconfined phase by their closed loop terms ( |5.11| ) 
with all possible windings while the magnetic correlators enter via perimeter 
contribution. As a result one has for the high-temperature phase (phase 2) 
(cf.pl). 

F{2) = l -eV z - {Nl - l)V 3 ^Q g - 7 -^N c V 3 T%Q q + O(N c ) (5.16) 

where Q q and fl g are perimeter terms for quarks and gluons respectively, the 
latter was estimated in ||157|| from the adjoint Polyakov line; in what follows 
we replace Q by one for simplicity. 
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Comparing (|5~T^| ) and (|5TTE|) , F(l) = F(2) at T = T c , one finds in the 
order 0(N C ), disregarding all meson and glueball contributions 



11 AT a °( p2 ) 

For standard value of G 2 = ^{F 2 } = 0.012 GeV 4 (note that for n/ = one 
should use larger value of G 2 |6£j) one has for SU(3) and different values of 
rif = 0, 2, 4 respectively T c = 240, 150, 134 MeV. This should be compared 
with lattice data ||152j| T l c at = 240, 146, 131 MeV. Agreement is quite good. 
Note that at large N c one has T c = 0(N®) i.e. the resulting value of T c doesn't 
depend on N c in this limit. Hadron contributions to T c are 0(N~ 2 ) and 
therefore suppressed if T c is below the Hagedorn temperature as it typically 
happens in string theory estimates 



Till now we disregarded all perturbative and nonperturbative corrections 
to F{2) except for magnetic condensate, the term ^eV 3 . If we disregard also 
this term, considering in this way only free gas of gluons and quarks for the 



phase 2, we come to the model, considered in [|159|| . The values of T c obtained 
in this way differ from ours not much - they are factor of 2 1//4 « 1.19 larger, 
but one immediately encounters problems with explanation of spatial string 
tension, screening masses etc., which are naturally are accounted for by the 
notion of magnetic confinement - nonzero values of magnetic correlators in 
the phase 2, including magnetic condensate term, l/2eV3. 

Notice, that NP corrections may contribute to fl g , fl q . Their phenomeno- 
logical necessity can be seen in the measured values of e — 3p, see the corre- 



sponding pictures in | |151| , |152|| . In case of fl g = Q q = 1 the difference e — 3p 
should be zero, and of course higher orders in N~ 1 (NP effects ) and higher 
orders in g (perturbative effects) contribute to it. In [ |157] ] the authors tried 
to estimate effect of nonzero (Q g — 1), which is 0(N 2 ), on the energy density 
and pressure. To this end the adjoint Polyakov line was used and the NP 
perimeter contribution was separated from it. 

5.3 Spacial W— loops 

In this section we derive area law for spatial W-loops, expressing spatial 
string tension in terms of CMC. 

In standard way one obtains the area law for large W-loops of size L, 
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L 3> T^ m ' [T^ 1 is the magnetic correlation length) 

(W(C)) 

spacial ~ GXp( ^s^mm 

where the spacial string tension is [|154| 



(5.18) 



2 

a s = 9 - J d 2 x((B n (x)$(x,0)B n (0)<£(0,x))) + 0((B 4 )) (5.19) 
and n is the component normal to the plane of the contour, while the last 



term in (|5.19| ) denotes contribution of the fourth and higher order cumulants. 
On general grounds one can write for the integrand in ( |5.19|) (see also section 
2): 

/ nD B \ 
{{Bi(x)$(x,0)Bj(0)$(0,x))) = <% \D B {x) + D B {x)+x 2 1 • 



dx 2 



gDf 

X^Xj— — „ . 

ox 



and only the term D B (x) enters in ( 5.19 ) 



2 

*' = jj d 2 xD B (x) + 0((B 4 )) 



(5.20) 



(5.21) 



similarly for the temporal W-loop in the plane (z4) one has the area law for 
T < T c with temporal string tension 



°e = ^J d 2 xD E (x) + 0((E i )) 



(5.22) 



For T = due to the 0(4) invariance CE and CM correlators coincide and 
oe = & S - The lattice measurements of De and Db |25], [2f| show that oe and 
os stay approximately constant in the whole interval < T < T c , as it is 
suggested by the phase transition mechanism described above. For T > T c 
in the deconfining phase CEC vanish, while CMC stays nonzero and change 
on the scale of the dilaton mass ~ 1 GeV, therefore one expects that a s 
stays intact till the onset of the dimensional reduction mechanism. This 
expectation is confirmed by the lattice simulation — a s stays constant up to 
T ~ 1.4 T c |160 . Lattice data ||160| , |161|| show an increase of a s at T rs 2 T c , 
for SU (2) which could imply the early onset of dimensional reduction. 
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It was shown how the Feynman-Schwinger representation (FSR) which 
proved to be very useful for T = 0, is modified for T > 0. In particular, 
all gluon and quark propagators are written as sums of path integrals with 
winding paths in the 4th direction. Using FSR the Hamiltonian for screening 
states of mesons and glueballs is derived, in |162| and screening masses and 
wave functions for lowest meson and glueball states are numerically com- 
puted. All this exemplifies NP background field theory and in particular 
MFC as a powerful instrument both for T < T c and for T > T c . 



6 Conclusions 

The material of previous chapters contains only a small part of results ob- 
tained with the method of FC. The aim of the present review is to explain 
foundations of the method and possible lines of development. 

In all fields where NP configurations are important the method of FC can 
be used to give a systematic description and calculation of effects. Here is 
the list of some results obtained heretofore with FC: 

• Proposed dominance of the lowest (Gaussian) FC gives a clue to the 
simple stochastic picture of confinement in agreement with lattice data. 

• The qq and gluon - gluon interaction is approximately proportional 
to the charge squared (quadratic Casimir operator C$) which is sup- 
ported by lattice data. 

• Heavy quark dynamics (masses, decay rates etc) calculated via FC are 
in good agreement with experiment. 

• The string profile (field distribution inside the string ) is well described 
by the lowest Gaussian correlator, for both types of probes (see chapter 
3). 

• The consistency of the model with low-energy theorems is established. 

• Lattice calculations yield FC with magnitude and correlation length 
T g which quantitatively correspond to the string tension and gluonic 
condensate known from experiment. The use of lattice measurements 
allows to put all the method of FC on the firm foundation, since using 
FC from lattice data enables one to predict hundreds of data: masses, 
widths, cross sections etc. 
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• It was predicted that at the deconfinement phase transition the corre- 
lator of CE fields D E vanishes, while that of CM fields D B stays intact. 
Lattice measurements have confirmed this and have shown the dra- 
matic decrease of De at T > T c , while Db is almost equal to its value 
at T = 0. From this fact it follows that De is a good order parameter 
for the deconfinement phase transition, and also that above T c one has 
the phase of "magnetic confinement". 



It is worth noting, that the fundamental input of the formalism - the 
Gaussian correlator - is not a freely adjustable function. Instead, one gets it 
from the lattice measurements and analytically calculates all the rest without 
any additional "fitting". It is one of the reasons why so much attention has 
been paid throughout this review to the question about Gaussian dominance, 
since in this case one basically does not need any other field-theoretical inputs 
(like higher correlators). 

One should also mention four topics which might have important influ- 
ence on the development of the method. An important recent development 
which has a direct bearing on FC is the calculation of the so called gluelump 
states on the lattice [119] and in the present method | 120|| . Actually gluelump 



Green's functions are generalizations of FC and the lowest gluelump states, 
1 and l + ~ refer directly to the FC of electric and magnetic fields, respec- 
tively. Therefore the masses of these gluelump states are simply T~ l for the 
corresponding correlators. 

The second line of development is the analytic approach with the aim 
of construction of the equations, defining the FC. In this way one should 
analytically calculate FC and establish the hierachy of FC corresponding to 
the Gaussian vacuum. Recently a set of integral equations has been derived 
in the large N c limit |70|| , and the first calculation of T g in terms of a was 
made possible, yielding results in good agreement with lattice data. More 
extended analysis is necessary in this direction including phase transition and 
dynamical quarks. 

As mentioned in the introduction, the model has also been successfully 
applied to high energy reactions with small momentum transfer. Here it was 
possible with the correlator discussed here to calculate the matrix elements 
of many high energy reactions involving hadrons and (virtual) photons. A 
discussion of the methods and the results can be found for instance in the 



reviews 



jql , |p7| and in and the literature quoted there. Also this is a wide 



field and will be treated separately. 
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Finally, there is another area of research important for applications and 
understanding - the interconnection of FC on one hand and OPE and QCD 
sum rules on another hand, especially from the point of view of OPE vio- 
lation. As it was argued by one of the authors [|133|| the problem of OPE 
violation is tightly connected to a vast and almost unexplored field of interfer- 
ence between purely perturbative and purely nonperturbative contributions. 
This direction of studies looks promising however difficult. The work on this 
set of questions |32], [33] in the example of e~e + annihilation raises many in- 
teresting problems, especially about ultraviolet renormalons and 1/Q 2 term, 
which together with recent lattice measurements of the condensate (FF) (see 
discussion and refs. in RB71) makes this field highly attractive and dynamic - 
hence still inappropriate for the review. 
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Figures 
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Fig.l. Df/a 4 = D + D 1 + z 2 dD l /dz 2 and D ± /a 4 = D + D 1 versus z. The 
lines correspond to the best fit to ( p.!4j ) (from |69| ) . 
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Fig. 3. The magnetic longitudinal correlator versus distance, for different 
values of T/T c (from |69||). 



73 




2 4 6 8 10 12 14 16 

r/a s 



Fig. 4. Static potentials for sources in different SU(3) representations (from 

OH)- 
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Figs. 6, 7 Difference of the (non-perturbative) squared electric field parallel 
to the quark (x 3 ) axis (—(E^(x 3 , r)) q q_ vac ) in GeV/Fm 3 for quark separation 
R w = 0.7 Fm (Fig.6, left) and R w = 3 Fm (Fig 7, right); T g = a = 0.35 Fm 
(from H). 
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